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Résumé. — Une photopile en circuit ouvert se charge progressivement lorsqu'on la refroidit : 
la différence de potentiel peut dépasser 2,5 volts à 100 °K. Le mécanisme de ce phénomène est 
étudié en détail ainsi que les processus de décharge obtenue de trois manières différentes ; par 
court-circuit, par éclairement et par réchauffement. On conclut enfin à l’existence de deux formes 
de la couche d’arrêt, stables l’une (C) à chaud (7 > 210 ©K), l’autre (F) à froid (7 < 210 ©K). 
Le passage d’une forme à une autre se ferait par migration d’électrons d’un niveau à un autre. 


Abstract. — An open-circuit photocell becomes progressively charged when cooled. At 100 OK, 
the potential difference may exceed 2.5 V. This process is discussed in detail, as is also that of 
discharge, which may be obtained in one of the three following: ways : by short-circuiting, by 
illumination and by heating. In brief, two stable forms of barrier layer are possible : one (C) obtai- 
at PT > 2100K, the other (F) at temperatures below that threshold. One surmises that the 
transition from one form to the other might be through migration of electrons from one level to 
the other. 


I. Apparition spontanée, aux basses températures, d’une différence de potentiel aux bornes de la 

photopile. — L'emploi d’un électromètre Curie à 
isolement élevé a permis de mettre en évidence 
l'apparition d’un nouveau phénomène. Lorsqu'on 
refroidit une photopile à l’obscurité il se crée à ses 
bornes une différence de potentiel qui croît avec 
l’abaissement de température et se stabilise aux 
environs de 1,3 volt vers 100 OK. 

La courbe représentant la variation de cette 
différence de potentiel en fonction de la tempé- 
rature est tracée au cours d’un refroidissement sur 
300 200 100 la figure 1. 

En régime statique on obtient les valeurs du 


F16G. 1. — Différence de potentiel et charge apparaissant 


aux bornes d’une photopile par refroidissement. tableau ci-dessous : 
TABLEAU 1 
TK 300 275 ‘ 9250 225 200 175 150 125 100 75 
V (volt) 0,000 0,000 0,000 0,000 0,050 0,140 0,210 0,900 1,300 1,300 
Ces mesures nécessitent un grand nombre de L'isolement doit être très soigné : nous avons dû 


précautions (classiques d’ailleurs) vu la valeur très éliminer même le polythène dont la résistance 
élevée atteinte par la résistance interne de la superficielle diminue beaucoup par temps humide, 
photopile aux basses températures. et qui tend à garder des charges adhérentes. 
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Nous n’avons utilisé que plexiglas et téflon. 
D'autre part pour éviter de recueillir sur la photo- 
pile des courants parasites, son circuit a été entiè- 
rement isolé et placé à l’intérieur d’un blindage 
jouant le rôle d’anneau de garde et mis à la masse. 
Un pôle de la photopile doit être relié à la masse : 
cette liaison est faite en un seul point, sur l’élec- 
tromètre même. Dans ces conditions nous avons 
obtenu des courbes parfaitement reproductibles à 
des mois d'intervalle. D’une cellule à l’autre les 
différences sont minimes. 

Un certein nombre de particularités caractérise 
ce nouveau phénomène. Si l’on maintient la cellule 
à la température T qui a créé l’apparition d’une 
tension V, celle-ci reste stable un moment puis 


V (voit) 


€ (heures) 


( 1 2 


F1G. 2. — Décroissance en fonction du temps 
de la différence de potentiel. 


F1G. 3. — Variation de la différence de potentiel 
par réchauffement consécutif à un refroidissement. 


décroît très lentement jusqu’à zéro. La figure n° 2 
montre un exemple typique de ce comportement. 
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Dans certaines expériences nous avons constaté 
une constante de temps de 2 heures 30 minutes ce 
qui correspond à une résistance d'isolement 
de 3 x 101% Q, soit 30 millions de mégohms. 

Si l’on réchauffe la cellule à partir de 90 °K la 
courbe V — f(T) fait apparaître une dissymétrie 
par rapport au refroidissement : nous en étudierons 
la cause plus loin. La figure 3 montre l’ensemble 
aller et retour 300 °K à 90 °K. 

Ces divers aspects du phénomène nous ont amené 
à penser qu'il s'agissait en l’occurrence d’une 
création par le refroidissement de paires électron- 
trou, c’est-à-dire en fait d’une véritable charge du 
condensateur constitué par la photophile, charge se 
maintenant un certain temps grâce à la valeur très 
élevée de la résistance de fuite. Mentionnons tout 
de suite que la différence de potentiel ainsi créée est 
de sens opposé aux différences de potentiel d’origine 
photoélectrique. ; 

Pour vérifier cette hypothèse nous avons réalisé 
le montage de la figure n° 4. 


F1G. 4. — Montage de l’électromètre pour mesurer la charge. 


Un condensateur variable étalon C peut être 
branché aux bornes de la photophile par le moyen 
d’un relais I monté spécialement sur plexiglas. La 
marche de l'expérience est la suivante : la cellule 
est refroidie à une température T' et l’électromètre 
indique une tension V,. Le relais I est alors 
enclenché et la tension V indiquée par l’électro- 
mètre varie en sens inverse de la capacité C. Il est 
possible de trouver une valeur C, de la capacité 
pour laquelle l’électromètre indique V,/2. On cons- 
tate alors entre V et C la relation suivante : 


(C, +'C) V = Ci Po soit Vol: 1 (CIC 


Le tableau ci-dessous reproduit des résultats de 
mesure : 


TABLEAU 2 


C 0 
LA 1,310 


100 
4,120 


200 
0,960 


300 
0,850 


400 
0,762 


De ces valeurs nous tirons les deux constantes : 


Co = 568pF et Vo — 1,310 volt. 

Voir la courbe de la figure n° 5. 

Si ces mesures sont effectuées pendant la durée 
où la tension en l’absence de C resterait constante, 
l'opération peut être recommencée plusieurs fois, 
C débranché et rebranché, ete. 


500 
0,693 


1100 
0,445 


600 
0,632 


700 
0,585 


800 
0,543 - 


900 
0,503 


1 000 
0,472 


Ces résultats prouvent que au cours de ces expé- 
riences On a partagé de diverses manières une 
charge électrique constante Q — C, V, entre une 
capacité fixe C, (photopile plus conducteurs plus 
électromètre) et une capacité variable C. 

L’abaissement de température a donc provoqué une 
libération d'électrons dans le sélénium et une accumu- 
lation de trous dans la couthe d'arrêt constituée par le 
contact sélénium-métal. 
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500 


1000 


Fra. 5. — Détermination graphique de C, et de V. 


La quantité d'électricité mise en jeu se calcule 
aisément à partir des résultats précédents : on 
trouve ; 

Q = Co 


en u. 6. M. : 
OSSI 10 PE SMILE 


Dans ce système la charge de l’électron étant 
représentée par le nombre 1,60 10—2°, on en déduit 
que le nombre de paires libérées est égal à : 


n = (7,5.10—11)/(1,6.10—2%) — 4,7 10°. 


La surface active de la photopile utilisée étant 
s — 0,636 cm?, le nombre de paires créées par 
centimètre carré est égal à : 


NES 7110 686 2 MON, 550107 


Sachant que la capacité de la cellule seule est 
de 300 picofarads nous en déduisons d’une part le 
potentiel réel à ses bornes.en l’absence d’appareil de 
mesure : 


V = V,(568/300) = 2,5 volts 
et l'énergie emmagasinée : 


E = (1/2) (@°?[C) —1,4.10—? erg par centimètre carré. 


IT. Décharge par court-cireuit. — La quantité 
d'électricité présente sur la photopile peut dispa- 
raître sous l’effet d’un court-circuit. 

Ainsi il est possible au cours d’un refroidis- 
sement, de court-circuiter à diverses reprises la 
cellule, les tronçons de courbe V — f(T) obtenus 
sont simplement décalés les uns par rapport aux 
autres sans modification ainsi que le montre la 
figure n° 6. 

La cellule déchargée par court-circuit reste dans 
l’état particulier qui a permis la création de charges. 
On s’en aperçoit à l’apparition d’un nouveau phé- 
nomène lors du réchauffement. De nouvelles 
charges de même sens font monter la tension aux 
bornes jusqu’à 0,5 volt environ, puis cette tension 
diminue, s’annule, s’inverse et s’annule enfin. La 
figure n° 7 montre ces variations : 

Cette courbe est bien définie et peut être repro- 
duite autant de fois que l’on veut à condition de 
respecter le cycle thermique. Toutefois la cellule 
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étant refroidie et déchargée peut rester dans le 
même état (à condition de rester à la même tempé- 
rature) aussi longtemps qu’on le veut. Nous n’avons 


200 NS 100 


Courts -circuits 


F1G. 6. — Refroidissement 
avec 2 courts-circuits intermédiaires. 


. 7. — Variation de la charge par réchauffement 
après court-circuit prolongé, à 90 °K. 


pas constaté de différence mesurable entre une 
courbe deréchauffement tracéeimmédiatement et la 
même courbe tracée (pour un autre cycle ther- 
mique évidemment) après avoir maintenu la cellule 
pendant 9 heures à 90 °K. 


III. Décharge par éclairement. — Un aspect 
intéressant de ce phénomène est la possibilité de 
décharger la cellule en l’éclairant: Nous avons pu 
mettre en évidence deux lois régissant cette 
décharge. 

19 La décharge est totale pour une énergie lumi- 
neuse incidente qui pour une longueur d’onde 
donnée ne dépend pas du facteur temps : seul 
compte le nombre de photons incidents. 

20 Le nombre de photons incidents nécessaire 
pour obtenir une décharge complète est fonction 
de la longueur d’onde. Connaissant la puissance 
lumineuse en watts et le temps on peut en déduire 
l'énergie en ergs nécessaire à chaque longueur 
d’onde pour décharger la cellule, On peut repré- 
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senter la variation en fonction de la longueur 
d’onde À soit de cette énergie, soit du nombre de 
photons qu’elle représente : ce sont les courbes de 
la figure n° 8. 
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Fi. 8. — Distribution spectrale de l’énergie et du nombre 


de photons nécessaires pour ramener la cellule dans son 
état «chaud». 


Il serait intéressant de voir également comment 
évolue avec la longueur d’onde, le nombre de 
paires électron-trou créées par effet photoélec- 
‘trique. Malheureusement ce décompte est presque 
impossible à faire car aux extrémités du spectre à la 
température de 100 ©K la sensibilité est extra- 
ordinairement faible et dans la partie jaune-rouge 
les phénomènes de resensibilisation font croître la 
sensibilité très rapidement, ce qui rend illusoire 
tout calcul. 

Une autre particularité de l’action de la lumière 
est que la courbe de décharge en fonction du temps 
a une allure qui dépend de la longueur d’onde. 
L'influence de la durée étant nulle il est possible 
d’étaler ces courbes à volonté pour mieux les 
comparer. La figure n° 9 montre les principaux 
aspects. 
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0,85 à 1,04 


FiG. 9. — Courbes de décharge à 90 0K 
en fonction de la longueur d’onde. 


Ces différences nous ont suggéré l’existence de 
deux phénomènes se masquant l’un l’autre partiel- 
lement. Nous avons cherché à les mettre en évi- 
dence de la manière suivante. 

Nous avons indiqué dans une publication anté- 
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rieure que certains rayonnements (d'efficacité 
maxima vers 0,75 x) avaient la propriété de resen- 
sibiliser une photopile préalablement refroidie. 
Nous avons cherché à voir si une cellule dont la 
charge était juste neutralisée par un éclairement de 
longueur d’onde À et d’énergie convenable était 
dans un état dépendant ou non de À. Le résultat 
est le suivant : l’état final de la photopile dépend 
de la longueur d’onde qui a été utilisée pour la 
décharger. En effet, la courbe de charge obtenue 


au réchauffement diffère selon les cas. La figure 


n° 10 en montre les aspects caractéristiques. 


0,4 à 0,6p 0,65 à 0,80 


0,85 à 1,0 


F1&. 140, — Courbes de réchauffement consécutives 
aux décharges précédentes. 


Il semblerait donc que, au moins à l’approxi- 
mation de ces mesures, le domaine spectral de 0,65 
à 0,80 à trouvé comme étant le plus actif pour 
ramener la cellule à l’état « chaud » soit d’étude 
relativement simple, l’énergie nécessaire à la 
décharge de la cellule étant dans ce cas juste sufli- 
sante pour resensibiliser la cellule. Nous avons 
d’ailleurs vérifié, en envoyant une énergie insuf- 
fisante pour décharger la cellule, que celle-ci n’était 
pas resensibilisée totalement. 

Par contre pour les deux extrémités du spectre 
les deux phénomènes semblent dissociés, l'énergie 
nécessaire pour décharger la cellule étant insuf- 
fisante pour supprimer tout phénomène ultérieur. 

Il faut remarquer cependant l'indépendance des 
deux phénomènes : décharge et retour à l’état 
« chaud ». En effet sur une cellule déchargée par 
court-circuit, il faut la même énergie pour la 
ramener à l’état «chaud » que si elle n’avait pas été 
déchargée. Par ailleurs nous verrons que certains 
cycles thermiques peuvent faire apparaître des 
charges de signe contraire qui disparaissent éga- 
lement sous l’action de la lumière. Il semble done 
difficile d’admettre le processus simple de neutra- 
lisation des paires électron-trou thermiques par des 
paires électron-trou photoélectriques comme suf- 
fisant pour interpréter le retour de la photopile à 
l’état «chaud ». 


IV. Décharge par réchauffement. — Nous avons 
vu l'effet du réchauffement sur la charge de la 
cellule. Mais sur une cellule refroidie, puis ramenée 
à l’état « chaud » par éclairement nous pouvons 
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effectuer le cycle suivant : réchauffer la cellule en 
court-circuit jusqu'à une température T7, la 
remettre en circuit ouvert et la refroidir. Comment 
varie la charge © obtenue en fonction de T ? La 
figure n° 11 le montre. 


AC (0 RE 


Q (e/cm°) 


.200 250 


Fr. 11. — Charge créée par refroidissement consécutif 
à un réchauffement à la température T portée en abscisse. 


Nous avons constaté que là encore le temps 
n'intervient pas. Une cellule maintenue plusieurs 
heures en court-circuit à 208 CK par exemple ne se 
charge aucunement par refroidissement. De même 
une cellule maintenue en court-circuit plusieurs 
heures à 230 ©K prendra par refroidissement la 
même charge que si elle n’y était restée qu’une 
minute. 

Il est intéressant de noter que la température 
de 210 °K représente le départ de la courbe de la 
figure n° 1 comme celui de la courbe de la 
figure n° 10. 


100 200 


Fi. 42. — Cycle de réchauffement et de refroidissement 
à partir de 90 °K pour une cellule court-circuitée à 90 °K. 


Reprenons l’expérience de décharge par court- 
circuit suivie de réchauffement décrite en IT (voir 
figure n° 7). La cellule se trouve alors déchargée à 
une température de 180 °K environ. Un refroidis- 
sement fait alors apparaître un nouveau type de 
courbe : figure n° 12. 

Le phénomène est plus lisible en portant en 


CRÉATION DE CHARGES ÉLECTRIQUES PAR REFROIDISSEMENT 


685 


abscisse alternativement réchauffement et refroi- 
dissement (fig. 13). 


100 180 2,100 


Fic. 143. — Cycle de réchauffement et de refroidissement 
à partir de 90 °K pour une cellule court-circuitée à 90 °K 
avec abscisses croissantes puis décroissantes. 


La cellule se trouve ramenée à 90 ©K et 
déchargée : ce cycle peut être reproduit indéfi- 
niment, semble-t-il. Nous avons une succession 
de 10 cycles enregistrés sur une durée de 4 heures 
sans modification apparente. 

Mais si dans la phase de réchauffement on 
dépasse 200 °K la courbe de refroidissement se 
modifie et se rapproche d’autant plus de la courbe 


n° À que la température de réchauffement a été 


plus élevée. 

La figure n° 14 reproduit ainsi les courbes de 
refroidissement obtenues après des températures 
de réchauffement T7 respectives de : 1909, 2100, 
2120, 300 °K. 


T'=300°K 


& 
ë 
ES T'=215°K 
D 
ot T'- 212°K 
T'=200°K 
"=190°K 
0 
200 150 100 T°K 
Fic. 44. — Courbes de refroidissement en fonction de la 


température de réchauffement à partir de 900, 


Encore une expérience sur cette cellule. Repre- 
nons la courbe de réchauffement de la figure n° 7. 
Après la branche ascendante traduisant l’appa- 
rition de charges d’une polarité déterminée, nous 
avons une branche descendante qui prouve l’annu- 
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lation de ces charges. La cellule étant toujours dans 
la forme F, cette annulation doit être due à l’appa- 
rition de charges de signe contraire. Pour le vérifier, 
au cours du réchauffement, arrêtons celui-ci 
à 120 0K et court-circuitons la cellule. En repre- 
nant le réchauffement nous voyons apparaître les 
charges contraire que nous avions prévues : 
figure 15. 


Q|(e/cm°) 


Court-circuit 


100 150 


F1c. 15. — Mise en évidence 
de l’apparition de charges de signe contraire. 


Les charges sont alors de même sens que celles 
produites par l'effet photoélectrique normal. Et 
pourtant là encore l’arrivée d’un nombre conve- 
nable de photons fait disparaître ces charges et 
ramène la cellule à l’état C. 


Deuxième partie 


Vue d’ensemble sur les phénomènes aux basses 
températures. — Nous pouvons considérer qu’une 
photopile au sélénium possède aux températures 
voisines de la température ambiante un état stable 
que nous appellerons «état chaud » ou C. La courbe 
de sensibilité spectrale est bien définie et ne pré- 
sente que de légères modifications en fonction de la 
température. 

Toutefois sa résistance interne (et sa force élec- 
tromotrice à vide) augmentent rapidement lorsque 
la température décroit selon la loi simple 
R — Riexp (AE JKT) mettant en évidence un 
niveau d’impuretés d'énergie AE — 0,20 électron- 
volt environ. à 

Lorsqu'on refroidit la cellule il apparaît progres- 
sivement à partir de 210 ©K un nouvel état que 
nous appellerons «état froid » ou F dont l’existence 
se manifeste de plusieurs manières : 

a) la sensibilité spectrale commence à baisser 
vers 210 OK ; 

b) la résistance À augmente toujours ; 

c) en circuit ouvert des charges apparaissent aux 
bornes de la cellule. 

Vers 100 0K la cellule s’est stabilisée dans la 
forme F qui n’est plus modifiée par un refroidis- 
sement plus poussé. Cet état est caractérisé par les 
phénomènes suivants : 

a) la sensibilité est pratiquement nulle, infé- 
rieure à 1/3 000 de ce qu’elle est à 300 K ; 
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b) la résistance s’est stabilisée à une valeur très 
élevée, ne dépendant pratiquement plus de la tem- 
pérature ; 

c) la chargée libérée atteint environ 101 élec- 
trons par centimètre carré à 100 °K et ne varie plus 
avec T°. 

Entre ces deux températures 100 0K et 210 0K 
existe toute une série d’états intermédiaires fonc- 
tions de la température et ne dépendant pas du 
temps de séjour à cette température. Les chan- 
gements d’état accompagnant les changements de 
température se produisent en un temps inférieur 
à celui que notre appareillage nous permet d’appré- 
cier, c’est-à-dire selon les cas, de quelques secondes 
à une minute. Nous allons voir que ces états ont 
des propriétés particulières. 

Les formes C et F ont des conditions et des 
domaines de stabilité nettement différents et les 
processus de passage de l’un à l’autre ne sont pas 
tous réversibles. 

Les trois croquis de la figure n° 16 le montrent : 

forme C : stable dans le domaine 7° > 210 °K 

forme F : stable dans le domaine 7 < 210 ©K. 


SA log R 


15 35 10 15 


8 5 10 15 ô 10 


F1. 16. — Courbes caractéristiques des états « chauds » 
 (G) et «froid » (EF): 

a) sensibilité spectrale ; 

b) log de la résistivité ; 

c) quantité d'électricité libérée. 


Aux basses températures la lumière transforme 
rapidement la forme F en forme C stable et dont 
les propriétés en fonction de 7° sont réversibles à 
toute température à condition que les cycles ther- 
miques étudiés ne traversent pas la température 
de 210 0K. 

Le refroidissement d’une température supérieure 
à 210 0K à une température inférieure à 210 °K 
transforme partiellement la forme C en forme F : la 
transformation est totale si la température est 
descendue à 100 0K environ. 

Si à aucun moment la cellule ne reçoit de lumière 
son état est parfaitement défini à chaque tempé- 
rature et les courbes F des. figures 164 et 16h 
peuvent être décrites indifféremment dans un sens 
ou dans l’autre. La courbe 16c n’est pas réversible 
car aux charges libérées pendant le refroidissement 
s'ajoutent les charges libérées pendant le réchauf- 
fement qui ne sont pas opposées (voir fig. 3). 
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En somme les apparitions et disparitions de 
charges électriques sur les faces de la couche d’arrêt 
ne sont autre chose que des témoins de transfor- 
mations se produisant au sein de la couche d'arrêt 
par effet thermique. 

Deux formes stables sont mises en évidence 
l’une C pour 7 > 2100K l’autre F pour T'< 210 0K. 
Le domaine 400 ©K < T < 210 °K est un domaine 
de transformations thermiquement réversibles. A 
une température quelconque l’action de la lumière 
transforme de manière irréversible la forme F en 
forme C qui reste stable à la température en 
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question ou au cours de cycles thermiques ne fran- 
chissant pas la température critique de 210 ©K. 

Cette transformation étant d’origine thermique, 
il est logique de supposer qu’elle est due au dépeu- 
plement d’un niveau au profit d’un autre qui se 
sature. 

Des photons d’énergie convenable (0,74 1 soit 
1,68 eV) sont susceptibles de faire réoccuper les 
places vacantes et de ramener la cellule dans son 
état initial. Certains niveaux seraient donc métas- 
tables. 

Manuscrit reçu le 18 avril 1958. 
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LARGEUR DE RAIES DANS LE SPECTRE DE VIBRATION-ROTATION 
POUR HCI ET CH,. VARIATION AVEC LA PRESSION 


Par Mme Mrrgize CAMEO, 


Laboratoire des Recherches Physiques de la Sorbonne, Paris. 


Résumé. —- Une méthode indirecte nous a permis la mesure des largeurs de raies dans le spectre 
moléculaire de vibration-rotation, avec un spectrographe Perkin-Elmer 12 C à prisme de LiF 
dont la faible dispersion ne conviendrait pas à des mesures directes. 

Après une vérification de la validité des formules de calcul et des mesures à pression variable, 


nous avons effectué : 


1) une étude de la variation, en fonction du nombre quantique de rotation, de la largeur de 
raies, pour la bande fondamentale de l’acide chlorhydrique. Cette variation s’avère en accord 


avec la théorie ; 


2) la mesure de la largeur des raies + 9 et + 10 de la bande v, du méthane. 


Abstract. — An indirect method for measuring the line widths,in the vibration-rotation mole- 
cular spectra, was applied to HCI and CH,, from measurements on a Perkin-Elmer spectrograph 


12 C with LiF prism. 


After establishing formulae for calculation, verification of their validity and a study of pres- 


sure broadening we have made : 


1) a study of line width dependance with rotational quantum number, for the fundamental 


band of HCI. 


The experimental curve obtained corresponds with theory ; 


2) measurements of line widths + 9 and + 10 of the CH, v; band. 


Introduction. — Il y a quelques années, la 
mesure de la largeur des raies de vibration-rotation 
de molécules simples à l’état gazeux se heurtait à de 
nombreuses difficultés expérimentales. 

La largeur réelle d’une raie représente une gran- 
deur difficilement mesurable directement sur un 
spectre et la valeur obtenue ainsi dépend, entre 
autres facteurs, de l’intensité de cette raie et de la 
largeur des fentes du spectrographe. Une mesure 
directe nécessite un système très dispersif permet- 
tant une bonne séparation des raies, une source 
très intense et l’utilisation de fentes très fines. 

Il paraît donc inutile maintenant de tenir 
compte, à cause du manque de dispersion des spec- 
trographes, des mesures directes effectuées avant 
1939 par des auteurs tels que Becker, Grasse, Chao, 
Margenau, Cornell, Watson, Herzberg, Spinks [1, 2, 
3, 4, 5, 6, 7, 8] sur des molécules gazeuses simples : 
HCI, HCN, NH: Plus récemment, Adel et 
Barker [9, 10, 11], grâce à un spectrographe très 
dispersif, ont pu obtenir directement la largeur des 
raies du gaz carbonique et de l’eau. 

Les méthodes indirectes anciennes ne paraissent 
pas non plus apporter beaucoup de précision à ce 
sujet. Par contre Penner et Weber [12], après sup- 
pression de la structure fine par addition d’un gaz 
inerte à pression assez élevée, mesurent la largeur 
totale de la bande, pour l’oxyde de carbone et 
Poxyde azotique, et en déduisent celle des raies. 
Cette méthode possède l'inconvénient de supposer 


égales les largeurs de toutes les raies, ce qui & priori 
n’est pas exact. 

La méthode indirecte de Mme Vincent [13] nous 
parait d’une utilisation plus générale et nous 
Pavons appliquée à la mesure des largeurs des raies 
de la bande fondamentale de l'acide chlorhydrique 
et de la bande v, du méthane, à l’état gazeux. 

Relativement aux corps que nous avons exa- 
minés, en laissant de côté les travaux anciens, nous 
trouvons des calculs théoriques de Lindholm [14] 
et de Foley[15] sur la largeur des raies de HCI, ainsi 
que le travail très important publié tout récem- 
ment par Benedict, Herman, Moore et Sil- 
verman [16] sur le même sujet (nous y reviendrons | 
plus loin). 

Les résultats sur la largeur des raies du méthane 
étant inexistants à notre connaissance, nous 
reprendrons l’étude commencée par Mme Vincent 
en 1954 [13]. 


PREMIÈRE PARTIE 


À. Méthode et mode opératoire. — Nous appel- 
lerons : 

À (en cm1) la demi-largeur de la raie à mi- 
hauteur, 

« (en em?) — /p l'intensité totale de la raie [p 
étant la pression, Z (en em? atm-1) l'intensité par 
unité de pression]. 

Nous définissons z une absorption totale, égale 


N° 8-9 


au produit de la largeur de fente dans le spectre a 
par l'absorption dans les conditions expérimen- 
tales 1 — T : 


= ai —T). 


Il existe une relation entre z, «, À, très complexe 
dans le cas général, qui se simplifie si la fente 
spectrale recouvre entièrement la raie et si l’absorp- 
tion n’est ni trop grande ni trop faible [13]. 


2? = rA pl (1) 


La technique expérimentale consistera à mesurer 
la transmission, 7, au maximum d’absorption des 
raies, la largeur des fentes dans le spectre aux 
diverses longueurs d’ondes, 
bante, {, et la pression, p. 

Tous les spectres ont été enregistrés sur spectro- 
graphe Perkin-Elmer modèle 12 C équipé d’un 
prisme de LiF, à quadruple passage. 

La méthode des lignes de base nous permet de 
déterminer T après enregistrement du spectre de la 
cuve vide et de la cuve pleine du gaz ; et selon le 
mode opératoire utilisé par Mme Vincent [13], nous 
calculons la largeur de fente dans le spectre pour la 
longueur d’onde de chacune des raies du gaz et 
pour chaque valeur de la fente utilisée. 


B) Études préliminaires. — Nous avons voulu 
tout d’abord effectuer, avec diverses cellules 
remplies d'acide chlorhydrique à la pression atmo- 
sphérique, une vérification de la formule (1). Elle 
consistera à étudier la variation de la valeur z en 
fonction de:la largeur f des fentes du spectrographe 
et à examiner la variation de z? lorsque nous utili- 
sons des épaisseurs absorbantes variables L. 


INFLUENCE DE f. — Pour obtenir une valeur 
correcte de z, donc de À, il faut que la largeur de 
fente dans le spectre soit inférieure à la demi- 
distance entre deux raies consécutives. 

Pour HCI nous avons utilisé quatre valeurs de la 
fente différentes, toutes inférieures à la limite infé- 
rieure que l’on calcule comme égale à 230 u, et les 
courbes z — /(f) tracées, par exemple pour les 
raies, + 2 et — 5, nous montrent que, si pour les 
valeurs faibles de f la variation de z est importante, 
elle diminue beaucoup pour des fentes plus larges 
(Ag. 1). 

En utilisant par la suite des fentes de 200 y et 
220 y, nous pouvons supposer la largeur de fente 
suffisamment importante pour couvrir la totalité 
de la raie et nous donner une mesure de l’absorption 
totale «vraie ». 

Des expériences similaires effectuées avec le 
méthane nous indiquent une fente utilisable com- 
prise entre 100 y et 130 y. 


a, l'épaisseur absor-: 
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INFLUENCE DE /. — Les courbes z? — f(l), tracées 
ici par exemple pour les raies intenses, — 2, et, 
— 5, s’avérent, aux erreurs d'expériences près, 


105 135 167 197 


F16. 1. — HCI. Variation de la quantité z 
en fonction de la largeur des fentes du spectrographe f. 


linéaires et passent par l’origine (fig. 2). Comme 
prévu, pour les raies très faibles, 9 et 10, où nous 
ne sommes plus dans des conditions suffisantes 
d’approximation, nous n’obtenons plus de droites. 


10,02 18,9 


Fi6. 2. — HCI. Variation de la quantité z° 
en fonction de l’épaisseur absorbante L. 


Dans cette première partie, nous avons défini une 
méthode et un mode opératoire, vérifié la validité 
d’une formule et étudié ses limites d'application. 

Il nous reste maintenant à aborder la mesure 
systématique de À et l'étude de sa variation en 
fonction de la pression pour l'acide chlorhydrique 
et le méthane. (Les valeurs de / n’étant pas obte- 
nues de la même manière pour les deux gaz, nous 
traiterons le problème séparément.) 
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DEUXIÈME PARTIE 


À) Acide chlorhydrique. — Nous rappelons les 
conditions opératoires : 
151,89; cm f = 200 up et 220 u. 
p compris entre 0,058 atm. et 1 atm. 


Étude de la bande fontamental de 3,46 u. 

Les mesures les plus récentes de l’intensité de 
cette bande, effectuées par Benedict, Herman, Moore 
et Silverman [16] donnent : Z, (intensité totale de 
la bande) — 130 em? atm-1. Cette valeur bien 
qu'inférieure à celle que donnent Penner et 
Weber [17] en 1953 est confirmée par des résultats 
obtenus par Coulon, Galatry, Oksengorn, Robin et 
Vodar [18] en 1954. 

La formule (11) donne l'intensité Z de chaque 
raie en fonction de /, et de J le nombre quantique 
de rotation 


TEUTE Rare 1)J (II) 


vo Q 
avec 
h? 
RER es 
” 8x? MAT 


(h, K,T ayant les significations habituelles et M 
étant le moment d’inertie). 

Q représente la fonction de partititon de rotation 

Fe 1)J m1 ee 
Q= D 7 +1) += 

Nous avons calculé les valeurs de xA aux diverses 
pressions et un premier examen des résultats nous 
montre une augmentation régulière en fonction de 
la pression. En traçant les courbes 7rA — j(p) 


His 


F1G. 3. — HCI. Variation de rA en fonction de la pression p. 


(fig. 3) nous voyons que cette variation est linéaire. 
Nous remarquons toutefois des points aberrants 
pour les deux pressions les plus élevées et pour les 
raies les plus fortes. Dans le cas d’une pression 
élevée, une raie intense occupe un grand domaine, 
que ne recouvre pas entièrement la largeur & de la 
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fente dans le spectre. Ce domaine se définit comme 
l'intervalle de fréquence dans lequel T (la trans- 
mission) se trouve inférieur à 1 — s(epetit)on a 


/Ip 1 pe 
d = 0 A (en cm): 


On obtient pour la raie + 3 et à la pression atmo- 
sphérique d — 10 em-1, alors qu’à cet endroit 
a = 5 cmt, La valeur exacte de xA, à la pression 
atmosphérique pour les raies les plus intenses, sera 
une valeur extrapolée sur les droites tA — f(p) 
tracées grâce à des mesures aux faibles pressions. 

Ces valeurs de xA ainsi obtenues ne tiennent pas 
compte de la présence, dans le gaz utilisé, des deux 
isotopes #CI et ?7CI qui donnent les composés H 35C] 
et H %7CI, Chaque raie de rotation, que nous avons 
mesurée se compose en réalité d’un doublet séparé 
par 2 cm1; et la largeur de fente utilisée ne 
permettait pas de les séparer. Or pour les raies 
intenses, nous avons vu que le domaine occupé par 
la raie peut atteindre 10 cm1; dans ce cas les 
deux raies empiètent totalement l’une sur l’autre 
et la largeur de chacune des composantes sera égale 
à la valeur calculée plus haut. 

Nous admettrons ce qui précède pour les raies 
— 5, —4, —3, —2, —1, +1, +2, +3, +4, 
+ 5, + 6. Par contre les raies- faibles — 10, 
— 9, — 8, +9, + 10 peuvent être supposées 
parfaitement séparées : dans ce cas la formule (III) 
donne À réel en fonction de la valeur calculée 


av” %1 Lo ne 
dy + 

Le rapport des intensités des deux isotopes &; /@» 
est environ 3. 


Acalc. = TE SA (I) 


À cm”! 


A réel = f(d) 


10 9 8270605 409 2010041278: 47 ENTRER 


Fic. 4. — HCI. Variation de la demi-largeur de raie 
en fonction du nombre quantique de rotation J. 


(Nous avons évidemment supposé égales les 
largeurs de chacune des composantes.) 
Nous pouvons dès lors tracer une courbe (fig. 4) 
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donnant les variation de la largeur réelle de chaque 
raie en fonction du nombre quantique de rotation 
et donner le tableau de ces diverses valeurs 
(tableau T). 


INTERPRÉTATION. — 1) Les courbes xA — f(p) 
passent par l’origine, ce qui correspond à une valeur 
négligeable de la largeur naturelle de la raie et de 
et de l'élargissement par effet Doppler par rapport 
aux autres causes d’élargissement. 


2) Les courbes tA —f(p) sont linéaires, ce qui 
nous amène à conclure à un élargissement uniquement 
dû aux chocs des molécules avec leurs voisines. 


3) Il existe une très grosse variation de la largeur 
des raies en fonction du nombre quantique de rota- 
tion, variation difficilement explicable si, au cours 
des chocs, toutes les molécules agissent de façon 
identique. Mais pour HCI, divers auteurs, 
Cornell [8] par voie expérimentale, Foley [15] et 
Lindholm [14] par leurs calculs et, plus récemment 
Benedict, Herman, Moore et Silverman [16], sont 
d’accord pour dire que tous les chocs ne sont pas 
identiques. 

Le nombre des molécules dans un état de rota- 
tion donné est proportionnel à la distribution ther- 
mique de Maxwell-Boltzmann. Un phénomène de 
résonance se produit lors du choc de deux molécules 
dans des états de rotation voisins et cause un élar- 
gissement supplémentaire. Ainsi plus les molécules 
seront nombreuses dans un état de rotation donné, 
plus elles auront de chances de rencontrer des 
molécules dans des états voisins. 

Lindholm [14] calcule ainsi l’élargissement par 
résonance lorsque J'— J = + 1 pour toutes les raies 
de la branche négative et de la branche positive. Il 
obtient un effet symétrique pour les raies des deux 
branches avec un maximum d’élargissement pour 
les raies 2, 3, 4, 5. Puis il fait le calcul pour 
J'— J —°, dans la branche négative et pour 
J'— J = ÿ? dans la branche positive et il obtient, 
dans ce cas, un élargissement supérieur pour la 
branche négative. Pour |[J' — J| > + 2, on peut 
négliger l’effet de résonance. Ces calculs expliquent 
l’allure de notre courbe À — f(J) en précisant que 
nous n'avons fait que mesurer la largeur totale due 
à l’addition des diverses causes d’élargissement. 


LARGEUR DE RAIES DANS LE SPECTRE DE VIBRATION-ROTATION 


691 


B) Méthane. — Soit ! — 1,89 cm, f — 104 u 
et 124 u et p compris entre 0,057 et 0,708 atm. 

Tout d’abord nous n’étudierons pas la variation 
de xA en fonction de p, mais celle de la quan- 
tité z?/pl qui, pour une même raie, se comporte en 
fonction de p, de la même manière que A et qui 
ne fait pas intervenir la valeur de / très contro- 
versée (fig. 5). 


0,50 


*/Plef(p) 


raie + 10 


? atm 
0,708 


G106 019% 0400 5% 


Fi. 5. — CH4. Variation de la quantité z?/pl 
en fonction de la pression p. 


Contrairement à ce que nous observions pour 
HCI, l'augmentation de cette quantité, d’abord 
très rapide pour les faibles pressions, diminue pour 
des pressions plus fortes. En traçant les 
courbes z? — f(p) nous voyons que la variation 
d’abord lente devient rapidement linéaire (fig. 6). 
Nous nous servirons de ces deux séries de courbes 
pour extrapoler une valeur de z? à la pression 
atmosphérique où, en raison de difficultés expéri- 
mentales, nous n’avons pu effectuer de mesures 
directes. 

Ensuite à l’aide des valeurs de Z mesurées pour 
chaque raie par Mme Vincent pendant notre 
travail (mesures non publiées), nous calculons les 
valeurs rA aux diverses pressions. 


INTERPRÉTATION. — Les résultats ainsi calculés 
nous paraissent surprenants : ainsi pour les raies 


TABLEAU 1 


HCI. VALEUR RÉELLE DE LA LARGEUR DES RAIES DE CHAÇUNE DES COMPOSANTES 
EN FONCTION DU NOMBRE QUANTIQUE DE ROTATION 


J — 10 9 8 7 6 5 4 3 2 ( 
A em—1 0,05 0,065 0,09 0,11 0,1% 0,16 0,17 0,17 0,17 0,15 
J ol 1 2 + 3 Se # + 5 RME E à 4 6 + 9° ‘+10 
A em é 0,15 0,14 0,14 0,14 0,135 0,11 0,10 0,07 0,05 0,0% 
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d'ordre J élevé, la valeur xA devient quatre fois 
supérieure à celle de raies plus étroites. 

L'étude effectuée par Boyd, Thompson et 
Williams [19] nous apporte l’explication de ce phé- 


p atm 
0,708 


0,06 0,194 0,400 0524 


F16. 6. — CH,. Variation de la quantité 7°? 
en fonction de la pression p. 


nomène par l'enregistrement d’un spectre du 
méthane à une très grande dispersion et par l’expli- 
cation qu'ils donnent de sa complexité : 

« Chaque raie peut se partager en une ou plu- 
sieurs composantes à cause d’interactions de 
Coriolis. On peut observer en outre la présence 
effective des raies d’une bande de vibration en 
principe interdite. De plus à la partie négative de la 
bande viennent se mêler les raies dues aux 
bandes 2y, et vo + va.» 

Il en résulte que nous avons probablement 
mesuré parfois la largeur d’un ensemble de raies, 
à cause de leur très important enchevêtrement. 

Après élimination des résultats obtenus : 1) pour 
les raies simples + 6, + 7, + 8, à cause de leur 
trop grande intensité, 2) pour les raies + 1 et + 2 
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à cause de la proximité de la branche Q, nous 
avons choisi deux raies simples + 9 (3 113,8 em=t} 
et + 10 (3 122,9 cm1). Leurs domaines dans le 
spectre n’empiètent pas sur celui de leurs voisines 
ce qui conduit aux largeurs 


NE, =" 000 eme 
NE =—10 0887 


à la pression atmosphérique 

Ces valeurs, faibles, surtout si l’on admet 
l'élargissement dû aux choc accrus par des 
actions intermoléculaires, dépendent de l’intensité 


_des raies. Welsh, Pashler et Dunn [20] qui donnent 


une valeur 1, (intensité totale de la bande) infé- 
rieure à celle que Mme Vincent a obtenue, ne 
paraissent pas tenir compte de la complexité de la 
bande. De plus si nous adoptions le chiffre des 
précédents auteurs À diminuerait encore. 

Pour expliquer cette faible valeur de la largeur de 
raie, nous avons effectué des mesures sur des 
spectres de mélanges de méthane et d’air, les 
pressions partielles du méthane étant les mêmes 
que précédemment et la pression totale corres- 
pondant toujours à la pression atmosphérique. Le 
tableau 2 donne les valeurs de z?/pl dans ce cas 
pour les raies + 9 et + 10 ainsi que celles qui 
avaient été précédemment trouvées. 

Au fur et à mesure que la proportion d’air croît 
la quantité z?/pl augmente. Donc, d’après nos 
expériences, l'élargissement causé par l'air est 
supérieur à celui que produisent des molécules de 
méthane, ce qui revient à faire intervenir, dans 
l'élargissement par chocs, le diamètre moléculaire. 
En effet, plus la taille d’un nombre donné de molé- 
cules reste petite, moins le nombre de chocs est 
important et moins l'élargissement est grand. Lei le 
diamètre des molécules de méthane, inférieur à 
celui des molécules d’azote, nous explique le phéno- 
mène d’accroissement de la valeur z?/pl lorsque la 
proportion d’air augmente, ainsi que la faible 
largeur des raies du méthane. 


Conclusion. — Nous pouvons comparer les résul- 
tats précédents avec ceux obtenus par Penner et 


TABLEAU 2 


CH4 VALEURS DE 2?/pl POUR DIFFÉRENTES PRESSIONS AVEC ET SANS ADJONCTION D’AIR 
POUR LES RAIES + 9 ET + 10. 


p (atm.) 0,708 0,525 0,400 
Raie + 9 0,91 0,96 0,98 

avec air 

sans air 0,81 0,73 0,65 
Raie + 10 

avec air 0,83 0,90 0,90 

sans air 0,73 0,69 0,60 


0,194 0,106 0,0703 0,057 
1,07 1,15 1,22 1,10 
0,54 0,40 0,36 0,20 
1,03 1,00 0,9% 0,98 
0,49 0,34 0,31 0,21 
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Weber [12] d’une part et Mme Vincent [13], d'autre 
part pour l’oxyde de carbone et l’oxyde azotique. 


Leurs résultats concordent et indiquent un élargis- 


sement indépendant du nombre quantique de rota- 
tion et linéaire en fonction de la pression. Ainsi 
Mme Vincent donne : 


0,038 cm—!1 
0,040 em—1. 


ÂAco = 
Axo = 


Dans ces deux cas particuliers et à la pression 
atmosphérique, la largeur de la raie est uniquement 
due au phénomène de chocs. L’élargissement iden- 
tique, qu'il soit dû à des molécules de la substance 
ou à celles de l’air, correspond bien à des diamètres 
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à peu près analogues pour ces deux gaz et égaux à 
ceux de l’azote. 

La valeur À ,,9 — 0,040 cm-! obtenue pour HCI 
nous montre que, pour des niveaux de rotation 
d'ordre plus élevé, lorsque le phénomène de réso- 
nance devient négligeable, la largeur de raie serait 
inférieure à celle de CO et NO ce qui correspond 
bien à un diamètre moléculaire plus faible de HCI. 

Ce travail a été exécuté sous la direction de 
M. Lecomte, Directeur de Recherches au C. N. R.S. 
et je tiens à le remercier par l'intérêt qu'il a pris à 
ce travail. Mes remerciements s'adressent aussi à 
Mne Vincent pour ses judicieux conseils. 


Manuscrit reçu le 18 avril 1958. 
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EMPLOI DE LA MASSE EFFECTIVE ET DE LA MASSE ÉLECTRONIQUE NORMALE 


DANS LES SEMI-CONDUCTEURS 


Par DANIEL CURIE, 


Laboratoire de Luminescence, Faculté des Sciences, Paris. 


Résumé. — J. Friedel a suggéré que pour décrire un niveau d’impureté,-une équation du type 
Wannier mais renfermant la masse ordinaire de l’électron libre peut donner dans certains cas une 
meilleure approximation que l’équation usuelle renfermant la masse effective des électrons de 
conductibilité. Ceci semble le cas pour les pièges à électrons dans certains cristaux luminescents, 
et en particulier dans les ZnS phosphorescents. 

Il faut pour cela que le développement en fonctions de Bloch de la fonction d’onde du niveau 
d’impureté ne renferme que des termes correspondant à des énergies suffisamment élevées dans la 
bande de conductibilité (électrons sensiblement libres) ; ceci peut être réalisé dans certains cas 
pour des niveaux profonds. Par contre dans le cas de niveaux peu profonds (donneurs dans Ge 
et Si), le développement ne renferme que des termes voisins du bas de la bande: on a alors l’équa- 
tion de Wannier avec la masse effective. 

Application (à défaut du réseau ZnS réel) au réseau de Kronig et Penney à une dimension. 


Abstract. — J. Friedel has suggested that in some cases an impurity level may be better des- 
cribed by a Wannier-type équation using the free electron mass than by the usual Wannier equa- 
tion using the effective mass. This is perhaps the case for electron traps in sulphide phosphors. 
Such a proceeding is valid for some deep levels, when, in the development of the wave function of 
the impurity level into series of Bloch functions, high-energy terms only appear : these terms may 
be associated with almost free electrons in the conduction band. On the other hand, shallow 
donors in Ge and Si are described by a development using only low-energy Bloch functions : thus 
the effective electron mass at the the bottom of the conduction band must be introduced into the 


Wannier equation. 


I. Introduction. — 1.1.— On sait que le calcul 
d’un niveau d’impureté par l’équation de Wannier 
n’est légitime que si ce niveau est très proche de la 
bande de conductibilité [1]. En désignant par Æ 
l'énergie du niveau créé par le potentiel pertur- 
bateur V; et par e — E — E, la profondeur de ce 
niveau au-dessous du bas Æ£, de la bande de conduc- 
tibilité (cf. fig. 2), cette profondeur est alors donnée 


par : 
ñ2 
— ON 
| 2m* 


(cas d’une bande isotrope non dégénérée ayant son 
minimum à k — 0); F(R,) est le coefficient de la 
fonction de Wannier w(r — R,) localisée autour du 
nœud R; du réseau dans le développement de la 
fonction d’onde d(r) décrivant le niveau d’impu- 
reté : 


+ Mn] F) = € FN) (1) 


dr) = Z F(R;) w(r — R)j). 


Si F(r) varie lentement on peut aussi bien écrire 
avec Luttinger et Kohn [2]: 


gr) © F(r). o(r) 


@o(r) étant la fonction de Bloch du bas de la bande 
de conductibilité. En ce cas le niveau d’impureté 
est donc décrit par la fonction F(r) modulée avec la 
période du réseau. 


The one-dimensional Kronig-Penney lattice model is studied as an example. 


Dans le cas d’un niveau situé profondément en 
dessous de la bande de conductibilité, il n’est plus 
permis de calculer l’énergie du niveau par l’équa- 
tion (1). J. Friedel a montré qu’on devrait en 
général obtenir une approximation plutôt meilleure 
en introduisant dans (1) la masse normale m de 
l’'électron libre, plutôt que la masse effective m* [3] : 

2? 


2m 


A + i(r)] Fr) = e Fr). (2) 


On a d’ailleurs l’idée intuitive qu'un électron 
situé sur un niveau d’impureté peu profond est 


presque libre dans une grande portion du réseau, —, 


d’où l'introduction de la masse effective corres- 
pondant au bas de la bande de conductibilité —, 
tandis que pour un électron situé sur un niveau 
plus profond ou plus localisé la masse effective des 
électrons de conductibilité n’a pas à intervenir (?). 

En général, l’approximation obtenue par la 
substitution pure et simple de m à m* dans l’équa- 
tion de Wannier{l’équation (2)]ne sera pas bonne ; 
elle pourra être meilleure si m* est déjà très voisine 
de m. k 


1.2. — Cependant le cas des pièges à électrons 
dans le sulfure de zinc [4] a donné un très bon 
accord entre les expériences (phosphorescence et 


(?) J. Friedel, S. Nikitine, Communications privées. 
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thermoluminescence) et cette façon de calculer la 
profondeur des niveaux par (2) et non par (1) : 

I] semble exister dans les ZnS phosphorescents : 

a) des pièges peu profonds, se laissant décrire 
par le modèle « en atome hydrogénoïde », c’est-à- 
dire dans lesquels l’électron capturé gravite autour 
d’une charge + e substitutionnelle ; la profondeur 
calculée pour ces pièges (qui s’identifie avec l’éner- 
gie d’activation thermique aux basses tempé- 
ratures [5]) a pour expression : 


c—RJKÈ  R—13,5 eV 


où X, est la constante diélectrique effective donnée 
par : 


1 d Das 1 
er ar ar 
Avec X — 8,3 constante diélectrique statique, 


K5 —5,07 constante diélectrique en haute fré- 
quence, on prévoit : 


Ecalculé — 0,28 eV 


ce qui coïncide exactement avec la valeur moyenne 
des résultats obtenus par les différents auteurs pour 
le piège responsable du pic de thermoluminescence 
vers 150 0K (cf. [4], Tableau I). 

Le même accord — surprenant et probablement 
fortuit dans une certaine mesure — est obtenu 
avec Cds : 


PIÈGES EN ATOME HYDROGÉNOIDE 


Ecalculé €observé 
Zn$ 0,28 eV 0,28 eV 
Cds 0,17 eV 0,17 eV 


b) des pièges plus profonds, associés au lumi- 
nogène, et Se laissant décrire par le modèle « en 
molécule HŸf », constitués par le piège précédent 
perturbé par une deuxième charge + e, due au 
centre luminogène voisin lorsqu'il est ionisé. La 
profondeur obtenue dépend de la distance r entre 
ces deux charges ; les calculs numériques peuvent 
être conduits jusqu’au bout, nous donnerons sim- 


plement la formule approchée : 
c— (RJKÈ) + (e?/Kr). 


Suivant les positions de voisinage entre ces deux 
charges + e, on obtient : 


PIÈGES EN MOLÉCULE H# Dans ZnS 
2 


€calceulé €observé 
0,815 eV 0,8 eV 
0:148/eN? (wurtzite seulement) 0,72 eV 
0,68 eV { 0,66 eV (souvent 
0,63 eV à non 
0,61 eV | 0,63 eV séparés) 
0,565 eV 0,57 eV 


Il existe vraisemblablement de tels pièges aussi 
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dans CdS, mais là on manque de données pour une 
comparaison avec l’expérience. 


1.3. MODÈLE HYDROGÉNOIDE DANS D'AUTRES 
CRISTAUX IONIQUES. — On observe un accord sem- 
blable : 


a) Oxyde de zinc. — Avec K — 12, K, — 3,88, 
on prévoit pour le modèle « en atome H » la pro- 
fondeur thermique 


£calculé = 0,26 eV. 


La thermoluminescence [6] donne deux pics, le 
plus important à 0,26 eV, un ie plus profond 
à 0,38 eV. 


b) Oxyde cuivreux. — Avec K —9, K, =4,0on 
a : 


NIVEAU ACCEPTEUR DANS Cu:O 


€calculé €observé 


0,32 eV 0,32 eV 


Il s’agit ici de l'énergie d’activation de la conduc- 
tibilité à l'obscurité [7]. 

Par contre : 

€) Pour PbS, avec K — 17,9 et K;. = 2,81, on 
prévoit un niveau hydrogénoïide à une profondeur 
de 0,30 eV, qui ne paraît pas avoir été observé. La 
largeur de la bande interdite étant d’environ 
0,34 eV, un tel niveau peut être noyé dans la 
queue de la bande d’absorption fondamentale. 


d) Dans les halogénures alcalins, le modèle hydro- 
génoide donne l’ordre de grandeur des énergies des 
centres F, mais les valeurs obtenues ont besoin de 
très fortes corrections, car l’électron capturé par 
une lacune d’ion négatif a, dans l’état fondamental, 
une grande probabilité de se trouver à l'intérieur 
de la lacune elle-même. 

Les valeurs 0,26 + 0,04 eV (dans ZnS irradié au 
Wood), 0,24 à 0,29 eV (dans ZnO) et 0,32 eV 
(dans Cuz0) ont été aussi obtenues par R. Frey- 
mann et ses collaborateurs comme énergies d’acti- 
vation de l’absorption dipolaire Debye [7]. Il faut 
cependant se méfier de telles coïncidences, comme 
Pa montré J. Le Fèvre (Diplôme d'Etudes Supé- 
rieures, 1957). 


[.4. DANS LE GERMANIUM ET LE SILICIUM, On à 
affaire à des niveaux beaucoup moins profonds et il 
est indispensable d'introduire la masse effective 
électronique dans l’équation de Wannier. Les 
valeurs expérimentales de l’énergie d’ionisation 
pour les impuretés P, As et Sb sont [8] : 


ÉNERGIES DES DONNEURS DANS Ge ET Si 


P As Sb 
Si 0,044 0,049 0,039 eV 
Ge 0,0120 0,0127 0,0096 eV 
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Avec l'emploi de la masse électronique normale, 
le modèle en atome hydrogénoïde donnerait, dans 
le silicium (X = K, = 11,7) : 


Ecalculé — 0,099 eV 
et dans le germanium (X = K, © 15,8) : 
€calculé — 0,054 eV. 


On trouve bien une profondeur de l’ordre des 
centièmes d’électron-volt, mais on ne peut guère 
parler d'accord avec l'expérience. 

Mais avec l'emploi de la masse effective m* (qui 
par suite de l’anisotropie est une moyenne conve- 
nable entre la masse longitudinale et la masse 
transversale), on doit multiplier les valeurs précé- 
dentes par le rapport m*/m et l’on trouve avec 
Kittel et Mitchell [9] : 

dans le silicium : 


* * à 
Mtransversale — 0,19 m Miransverante 0, J7077 
Ecaleulé — 0,0298 eV 

dans le germanium : 
* *k > 
Mtransversale — 0,082 m Mtraneverentel=- 0 7 


Ecalculé — 0,00905 eV. 


L'accord avec l'expérience se trouve grandement 
amélioré ; il reste encore de notables corrections, 
attribuées aux effets de cœur (?). Ces corrections 
sont plus fortes pour le silicium que pour le ger- 
manium, le rayon des orbites y étant moins 
étendu : 


ag 20 À AGe CES AU 


Dans ZnS on s'attend à des effets de cœur beau- 
coup plus considérables ; on a en effet : 


a = Ke @o ay = À? me? 


d’où a = 3,7 À : dans la wurtzite, l'orbite hydro- 
génoïde centrée sur une lacune soufre renferme en 
son intérieur le zinc premier voisin, à la distance 
d —2,35 À, mais passe sensiblement sur le zine 
suivant, qui se trouve à la distance 5 d/3. 

Dans NaCI, l'orbite hydrogénoïde centrée sur une 
lacune CI- aurait un rayon 4 © 2 À, tandis que le 
rayon de la lacune est 2,6 À environ : l'orbite pas- 


(2) Il faut raccorder le potentiel coulombien, valable à 
grande distance de l’impureté, avec le potentiel de Hartree 
valable dans la cellule où se trouve l’impureté : on abaisse 
alors l’état fondamental (c’est-à-dire qu’on augmente |<| 
d’une manière qui dépend de la nature de l’impureté [8]. 

Le problème envisagé dans ce travail est différent : il 
s’agit de savoir par quelle équation l’on devra remplacer 
l'équation de Wannier avec la masse effective (1) — 
valable pour un niveau peu profond et peu localisé — 
dans le cas d’un niveau plus profond, mais toujours peu 
localisé. Dans le cas d’un niveau plus fortement localisé, 
l'équation du type Wannier, soit (2), soit (1), ne fournira 
bien entendu qu’une première approximation sur laquelle 
il conviendra d’effectuer les corrections d’effet de cœur si 
l’on veut raffiner le modèle. 
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serait dans la lacune et cette fois le modèle hydro- 
génoïde n’est plus valable. 

Ainsi l’on ne serait attendu à un accord avec 
l'expérience bien plus médiocre pour (2) dans le 
cas de ZnS, CdS, etc. que pour (1) dans le cas de Si 
et Ge. Or c’est le contraire qui s’est produit. Peut- 
être nos approximations ont une validité meilleure 
qu’on le penserait à priori ; mais peut-être aussi ce 
bel accord risque de disparaître dans une théorie 
plus élaborée. 

L’étude du mécanisme de la phosphorescence, en 
particulier, montre que le « bas de la bande de 
conductibilité » est une notion fort mal définie dans 
un cristal réel ; il en sera donc de même de l’énergie 
d’ionisation d’un niveau d’impureté, surtout 
lorsque celui-ci est peu profond : nous pouvons 
concevoir que l’équation de Wannier (1), rigou- 
reuse pour un niveau peu profond dans le cristal 
idéalisé qui fait l’objet de nos calculs, se trouve 
donner dans le cristal réel de moins bons résultats 
que l’équation approchée (2) pour un niveau plus 
profond. 


IT. Conditions générales pour l’emploi de la 
masse normale au lieu de la masse effective. — 
I1.1.— Un niveau d’impureté peut être décrit par 
un développement selon les fonctions de 
Bloch o,(K, r), comprenant en toute rigueur tous 
les nombres d’onde K et toutes les bandes n : 


Yr) — 2 cu(k) qnlksr). (3) 


Nous ne pouvons arriver à une équation du type 
Wannier, soit (1), soit (2), que si nous avons affaire 
à un niveau localisé encore suffisamment proche 
de la bande de conductibilité, par exemple, pour 
que le développement (3) ne renferme pratiquement 
que les w, de cette bande. Cette hypothèse faite, 
on obtient l’équation de Wannier sous la forme : 


[ent iv) + Vi] Fr) = e Fr) (4) 


&,(K) est l’énergie de l’électron dans la bande de 
conductibilité, e la profondeur du niveau au- 
dessous de cette bande et Æ(r) la transformée de 
Fourier de c,(k). 

Lorsque le développement (3) ne renferme que 
des fonctions de Bloch du bas de la bande, on peut 
écrire : 

RE : ñ? 
Se A ha 7 2m* 


(bande isotrope à minimum pour # — 0) et l’on 
retombe sur l'équation usuelle (1). 

Par contre lorsque V; est assez profond on devra 
utiliser dans (3) des fonctions de Bloch qui corres- 
pondent à une énergie élevée dans la bande ; la 
condition pour la validité de (2) est que pour ces 
énergies l’on puisse écrire : 

en) = fiv) = — 5% A 


Dr 
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c’est-à-dire que l'énergie cinétique e,(k) dans la 
bande de conductibilité soit donnée par la même 
expression que pour les électrons libres, Cette con- 
dition ne saurait être réalisée que d’une manière 
approchée ; il peut suffire qu’elle le soit au voisi- 
nage du moment moyen k,. 

La validité de l’équation (2) nécessite plusieurs 
conditions simultanées : il faut que le niveau 
d’'impureté soit assez profond pour que k, soit 
élevé dans la bande, mais pas trop localisé cepen- 
dant, de manière que le développement (3) ne ren- 
ferme pas une contribution appréciable des ec, dont 
le K est voisin de zéro, ou d’une autre bande. 


11.2. EXPRESSION DU MOMENT MOYEN #,,. — On 
traitera |c,(k)[? comme la distribution de proba- 
bilité des diverses valeurs de k£. Considérons un 
niveau d’impureté du modèle hydrogénoïde : 


F(r) === exp (— r Ja) 


T 4° 
la transformée de Fourier de F(r) est : 


(2a)3/2 
SRE TE 


la distribution Z(k) des moments # est ainsi : 


32 %ak2 
2 Rs $ 
T(k) = 4x k?|c(k)| = +arh 


E(km)l en 6,78 | EI 


Bande ; 
de conductibilité 


Niveau d'impureté 


Fic. 14. — Énergie e d’un niveau localisé (hydrogénoïde) 
et énergie moyenne e(k») des fonctions de Bloch décri- 
vant ce niveau. Lorsque le niveau est peu profond, ces 
deux valeurs sont proches du bas de la bande de conduc- 
tibilité et l’on obtient l'équation de Wannier avec la 
masse effective. Si par contre le niveau est profond, 
e(kn) est elle-même de l’ordre de grandeur de e (cf. Cours 
de P. Aigrain, 1956-1957). 


Au moment moyen : 


8 


Pere 
d 3 ra 


correspond l'énergie emétique 


/ 
SE Le] — 0,218 [el 
QT 


(kw) = 


proportionnelle à la profondeur e du niveau. 
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On peut dresser le tableau suivant : 


c eV eff) eV 
Atome d'hydrogène dans le vide 19359 9,7 
Modèle H dans Ge 0,01 0,01 
Modèle I dans Zn$S 0,20 


0,28 


Admettre la validité de (2) dans ZnS et de (1) 
dans Ge revient à admettre qu'à une hauteur 
de 0,20 eV dans la bande, les électrons de conducti- 
bilité ont un comportement déjà beaucoup plus voisin 
de celui des électrons libres qu’à 0,01 eV : hypothèse 
plausible, mais on ne saurait affirmer qu’elle est 
réalisée qu’en présence d’un traitement théorique 
rigoureux du Zns, et alors les traitements approchés 
actuels du type Wannier seraient sans objet. 


II.3. ESTIMATION DE LA CONTRIBUTION DES 
DIVERSES BANDES DANS LE DÉVELOPPEMENT [3]. — 
Kittel et Mitchell en avaient déjà donné une pre- 
mière évaluation [9] ; une formule plus récente de 
Kohn [8]est : 

C E sir 


n 
RARE 
Cm . gap & 


ñn, désignant la bande de conductibilité et r une 
autre bande, par exemple la bande de valence. 
Comparons à cet égard Ge et ZnS : 


Gê Zn$S 
e profondeur du niveau 0,01 eV 0,28 eV 
Ezap bande interdite 0,7 eV 3,7 eV 
ar constante du réseau 5,62 À 511 À 
a rayon de l’orbite 5 NAS 3,7. À 
R rapport des amplitudes 0,0018 0,11 
[R|? rapport des probabilités 10% 10e 


L’admission d’états de la bande de valence dans 
le développement (3), bien que considérablement 
plus forte dans ZnS que dans Ge, y reste cependant 
une faible perturbation. 


III. Essais de justification pour l’introduction 
de la masse normale et difficultés rencontrées. — 
III.1. FORMULATIONS DIVERSES UTILISÉES POUR 
DÉCRIRE LE NIVEAU LOCALISÉ. — Nous commen- 
cerons par utiliser le développement en fonctions de 
Bloch (méthode d’Adams [1]) : 


Ur) = À en(k) qn(K, r). (3) 
Vu la relation entre les fonctions de Bloch On et 
de Wannier w, : 


En (K, r) = D eik.R; PA (T — R;) 
J 


il revient au même de poser : 


br) = Z F(R;) wifr — R;) 


698 
en définissant F(R;) par : 
D Cn(K) ER} 


n,k 
L’interpolation de F(R;) par une fonction #(r) 
qui est la transformée de Fourier de c,(k) ne pré- 
sente pas de difficultés si les sommations à effectuer 
s'étendent en fait à une seule zone (n — no). 
Si F(r) varie lentement nous pouvons sortir F(R;) 
de la somme précédente : 


F(R) — 


gr) = F(r) Zwnlr — R;) 


et comme en désignant par œ,(r) la fonction de 
Bloch du bas de la bande de conductibilité : 


Dolr) = Pnl0; Tr) — D es R;) 


nous obtenons : 
d(r) = F(r) polr). (5) 


Le niveau d’impureté se trouve alors défini par la 
fonction F(r) modulée avec la période du réseau 
(cf. Luttinger et Kohn, [2], [8]). Mais l’équivalence 
avec (3) n’est qu’approchée. 

J. Friedel utilise souvent le changement de 
variable (5), définissant par là la fonction F(r) [3, 
10]. Nous utiliserons successivement les deux for- 
mulations. 


111.2. MÉTHODE DU DÉVELOPPEMENT EN FONC- 


TIONS DE BLOCH. — Les +, sont les fonctions 
propres de l’hamiltonien du cristal parfait : 
H5 n = En(K) ?n (6) 


On(K,r) = walk, r) eik.r 


tandis que Ÿ est fonction propre de l’hamiltonien 
du cristal perturbé par le potentiel V; dû au 
défaut : 

(Hy+V)9—=EY (2) 


(les énergies seront désignées par Æ tant que 
l’origine des énergies n’est pas précisée, et par € 
lorsque le bas de la bande de conductibilité est pris 
comme zéro d'énergie). 

Nous serons amenés à postuler que la partie 
périodique u,(K, r) d’une fonction de Bloch corres- 
pondant à une énergie suffisamment élevée dans la 
bande de conductibilité est très peu modulée dans 
l’espace : 


on(K, r) : : eik.r pour k assez grand. 
L'exemple du réseau de Kronig et Penney nous 
montrera que cette hypothèse est assez plausible. 
Dans le domaine d’énergies où cette approxi- 
mation est valable, les électrons de conductibilité 
se comportent comme des électrons libres en ce 
sens que leur masse effective est égale à la masse 
normale de l’électron. Mais ceci ne suffit pas à 
assurer la validité de l’équation (2), il faut en plus 
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que le zéro d’énergies soit le bas de la bande de 
conductibilité : ceci n’a pas lieu dans le réseau de 
Kronig et Penney ; mais celui-ci est assez éloigné 
des réseaux cristallins réels. 

Soit VA le potentiel périodique du réseau : 


a? 
DS TA EPA 
Ho 2m pue 
L'énergie correspondant à une fonction de Bloch 
en er (comme V; n’est pas constant, ce ne peut 
être évidemment qu’une solution approchée) est : 


ES 2 p2 
Elk) re pds + [va dr. 
Le bas de la bande de conductibilité (4 — 0) sera 
le zéro d’énergies si 


Î Vr dr: 


En ce cas seulement, la profondeur € calculée 
par (2) sera à porter à partir du bas de la bande de 
conductibilité. Cette condition pourra être satis- 
faite, en particulier pour certaines propriétés de 
symétrie du potentiel Vx. 

Substituant le développement (3) dans l’équation 
d’onde (7), et tenant compte de (6) ainsi que l’ortho- 
normalité des fonctions de Bloch, on a l’équation 
bien connue (voir par ex. Kohn [8|) : 

2 (n, k|Viin’, K°) cv(k') = 0 (8) 


n°",k 


[En(k) — E] cn(k) + 
on à posé : 
(a, KPiin', &°) = fon (kr) Fi (kr) dr. 
Cette équation est rigoureuse. 


Si V; est constant, tous les c,(k) sont nuls sauf 
celui qui correspond au nombre d’onde K dans la 
nieme bande donné par : 


E,(K) E + V; = 0. (9) 


On a alors développé Ÿ avec une seule fonction 
de Bloch : 


Ÿ = Cn(K) Pn(K, Tr). 


Bien entendu si V; est vraiment constant, l’élec- 
tron n’est plus localisé ; en ce cas, on a affaire à un 
nouveau cristal parfait dont les bandes d'énergie 
sont décalées de V, : d’où ce résultat. 

Appelons Kk,, la valeur de K solution de l’équa- 
tion (9) : c,(k) a la forme d’une fonction 3(k —K,,) 
et sa transformée de Fourier F(r) est d'amplitude 
constante : 

F(r) : : eik 


in 2, 


On sait que F(r) est solution de l’équation de 
Wannier sous la forme (4). 

Cet exemple montre qu’un niveau d’impureté 
peut être à la fois profond (grand moment moyen k, 
si V; est grand) et peu localisé (faible étendue de la 
distribution c,(Æ) autour de k,) ; ces deux condi- 
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tions simultanées, énoncées au $ IL.1, sont ainsi 
parfaitement compatibles. 


Laissons maintenant V; quelconque, mais faisons 
sur les +, intervenant dans le développement (3) 
l'hypothèse u,(k,r) = const. : 


(a, KP n°, ke) = [ri ex ar, 


Bande de conductibilité 


Niveau 
localisé 


Potentiel perturbateur 


Fire. 2. — Lorsque le potentiel perturbateur V; est cons- 
tant, la fonction d’onde décrivant l’électron localisé fait 
intervenir la seule fonction de Bloch d’énergie e(k,,) dans 
la bande de conductibilité [3]. Si le niveau est profond, 
km sera grand ; si le niveau est peu localisé, il n’inter- 
viendra qu’une bande étroite d’énergies voisines de e(k,;). 


Introduisons là transformée de Fourier U de 
V; [21]: 
1 : 
— Jet V(r) dr 


nous reconnaissons l’équation de Wannier dans 
l’espace des moments : 


U(K) 


[Eutk) — E] cnlk) + er 2 U(k — k') c\(k) — 0 


: d’où l’on remonte immédiatement à (4) en posant. 


= 2 cn(k) eiker. 
Nous avons ainsi justifié (4) sans rencontrer les 
difficultés que produit, dans le cas général, l’intro- 
duction de la modulation ux (Kk, r) u,(k’, r) dans 
l'élément de matrice (| V;|) [2, 8]. Le mieux serait 
alors une connaissance complète de Æ,(k) pour 
toutes les énergies dans le cristal parfait, mais faute 
de cela nous pourrons appliquer l'équation (2) si : 


AR 


——_— 10 
2m (10) 


En(K) = 
au voisinage du momen moyen k,. 

Nous avons examiné les formes possibles pour la 
courbe e,(k) dans [11]: la difficulté n’est pas 
que £,(k) se raccorde avec la parabole des électrons 
libres, mais que cette parabole ait pour sommet le 
bas A de la bande de conductibilité (fig. 3). Si on 


MASSE EFFECTIVE ET MASSE ÉLECTRONIQUE NORMALE 


699 
admet pour le spectre des énergies Z,(k) l'allure 
classique figurant dans Seitz [12], ce raccordement 


serait plus aisé si la masse effective m* au bas de la 
bande de conductibilité est supérieure à la masse 


n (À) 
JÉ-RU4E (A) 


En (Am 
E(4) 


0 


ki Rs 


Fic. 3. — Courbe d’allure classique E,{k), 
(en trait plein) [12]. 

Les courbes en tirets représentent les paraboles des 
électrons libres, d’origines respectives O et le bas À de la 
bande de conductibilité. Elles se coupent en K. 

Sile moment k, est voisin de l’abscisse de K, l’approxi- 
mation de Wannier avec la masse normale est valable 
pour le niveau situé au-dessous de À comme pour le 
niveau situé au-dessous de ©. 

Pour la courbe en trait plein la courbure en A est 
faible et la masse effective au bas de la bande de conduc- 
tibilité plus forte que celle de l’électron libre; la 
courbe Æ,(K) passe par le point K et la méthode marche 
sous la condition précédente. 

En pointillé la courbe Z,(k) si m* < m ; elle passe ne 
pas par K et la méthode ne marche pas. Cette conclusion 
n’est valable qu’à une dimension. 

Si l’on accepte pour la courbe Æ,(k) une allure plus 
complexe que celle figurant sur ce Schéma (ce qui paraît 
bien être le cas en réalité), cette courbe pourrait coïncider 
avec la parabole des électrons libres issue de A sur un 
intervalle beaucoup plus grand, et la validité de la 
méthode serait générale pour un grand intervalle de 
valeurs de k». 


d’après Seitz 


normale M que si l’on avait m* << m (dans ZnS, 
l'absence d’effet d'extinction du champ magnétique 
sur l’électroluminescence suggère effectivement une 
masse effective élevée) (°). 


() Dans un cristal ionique tel que Zn$, la masse effec- 
tive est certainement beaucoup plus élevée que dans le 
Germanium. 

On a la règle de sommation : 


il 1 i 2 > [(n|pælno) 
mm m 


F7 En RCE ce 


où (ñ|pzlno) est l’élément de matrice de l'opérateur py 
entre les états du bas de la bande de conductibilité n, et 
de toute autre bande n. 

Dans Zn$ la bande de conductibilité est attribuée à des 
niveaux des ions zinc et la bande de valence à des niveaux 
des ions soufre ; ceci n’est qu’une approximation, valable 
dans un modèle ionique du sulfure, mais certainement 
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Il est possible que l'égalité (10) soit réalisée 
d’une manière beaucoup plus générale que le 
suggère le schéma de la figure 3. 

Si la masse effective est égale à la masse nor- 
male, m* — m, (10) sera automatiquement vérifiée 
tant que #,, n’atteint pas la limite supérieure de la 
bande. Si m* -£ m, le raccordement entre £,(k) et 
la parabole (10) aura lieu sur un intervalle de 
valeurs de £ d’autant plus étendu que m* sera 
voisine de "1. 


IIL.3. COMPARAISON AVEC LA MÉTHODE DE RAC- 
CORDEMENT INTERBANDES (J. des Cloizeaux et 
P. André [13]). 

Cette méthode repose sur l’équation générale de 

: Wannier : 


Eat—iv) + Vi] F(r) = EF(r). (4) 


Elle utilise des renseignements sur le compor- 
tement de Æ,(k) dans la bande interdite, fournis par 
prolongement analytique de ce qu’on sait expéri- 
mentalement sur la nature de cette fonction dans 
les bandes permises. 

Dans les bandes permises le nombre d’onde k 
est réel ; on connaît expérimentalement (dans Ge, 
pas dans ZnS$ !) les masses effectives m% dans la 
bande de conductibilité et m* dans la bande de 
valence. On peut donc poser : 

— au voisinage du bas de la bande de conduc- 
tibilité F4 : 

E = Eolk) = Enc +R 
2Me 

— au voisinage du sommet de la bande de 
valence E,, : 


E — Ey(k) ns Esy 


2? 
% 
2Mv 


k?. 


Dans la bande interdite À — 1 8 est imaginaire 
et Æ(k) se prolonge analytiquement suivant une 
fonction g(B). Au voisinage de E4e, On a : 


el n° 2 
g(B) = Ene OS 
2mMe 


et au voisinage de Æ4, : 


RUN 
g(B) = Enr + > É* 


2my 


exacte grosso modo. Dans ces conditions les fonctions 
d’onde du bas de la bande de conductibilité et du sommet 
de la bande de valence sont respectivement surtout loca- 
lisées, au voisinage des ions zinc pour la première, des ions 
soufre pour la seconde : les transitions entre les deux sont 
peu probables et contribuent peu dans la Y. (J. Friedel, 
Communication privée.) Le 

La masse effective m* pourrait être assez voisine de la 
masse normale dans la bande de conductibilité de ZnS, se 
rapprochant de celle du zinc métal. 

Si la contribution à la X des bandes plus élevées que la 

n 

bande de conductibilité (£, > £,,) l’'emportait sur celle 
de la bande de valence et des bandes inférieures (£y < Eh), 
m* pourrait même être supérieure à 7. 
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Si la bande interdite n’est pas trop large, comme 
c’est le cas pour le germanium, on peut déterminer 
par interpolation g(B) avec une précision très sufli- 


Fic. 4. — Raccordement interbandes [13]. 


En trait plein, la « courbe en forme d’œuf » g(f), pro- 
longement analytique du spectre des énergies Æ,(k) 
dans la bande interdite. 

En tirets, la parabole des électrons libres issue du 
bas À de la bande de conductibilité. 

Le cas de figure correspond à une masse effective mé 
au bas de la bande de conductibilité supérieure à la 
masse normale : la courbe g(B) et la parabole des élec- 
trons libres se coupent nécessairement en un point K 
Si kn est voisin de l’abscisse de K la méthode marche 
(application de l’équation de Wannier avec la masse 
normale). 

La « courbe en forme d’œuf » peut coïncider avec la 
parabole des électrons libres sur un plus grand inter- 
valle de valeurs de k, et cela même si mé < m, mais 
pour cela ilfaudrait accepter des formes plus compliquées 
pour g(f). 


sante (courbe en forme d'œuf, voir figure 4). 
Comme 
Eav)= EMI 
l’équation de Wannier s’écrira : 
(e— v) + Vi] F(r) = EF(r). 


Si la bande interdite est large, comme c’est le cas 
pour ZnS, il est à craindre que la donnée de m£ et 
mx ne suffise plus à déterminer g(8) ; et en parti- 
culier il peut arriver que g() soit voisin dans un 
certain intervalle de la courbe : 


> 

Cu 2 
D] 
2m 


g() = Eae ss (11) 
auquel cas l’on retomberait sur l'équation (2). 
Ainsi les deux méthodes sont complémentaires, 
mais il arrivera souvent qu'aucune des deux ne 
donnera mieux qu’un simple ordre de grandeur. 
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111.4. DESCRIPTION DU NIVEAU LOCALISÉ PAR 
UNE FONCTION DE BLOCH MODULÉE [2, 3, 10]. — 
Pour un niveau peu profond, on utilise (5), intro- 
duisant ainsi une fonction lentement variable Æ(r) 
modulée par la fonction de Bloch +, décrivant le 
bas de la bande de conductibilité. Pour un niveau 
plus profond, il est préférable d'introduire la modu- 
lation produite par la fonction de Bloch 


ok, Tr) = 1u(k;, r) eik,-r 


correspondant au moment moyen 4. 
Ainsi la fonction d’onde du niveau d’impureté 
sera posée : 


gr) = F(r) u(km, r). (12) 

On a modulé FÆ(r) par l’onde stationnaire u(ky, r) 
et non par l’onde progressive exp (1k,.r). 

Pour simplifier, nous allons traiter le problème à 
une dimension seulement ; à 3 dimensions, on peut 
être amené à introduire toutes les ondes de Bloch 
dont le moment k,,, dans un cristal isotrope, aurait 
la même valeur mais des directions différentes : 
posons : 


kr = kyt 
(t vecteur unitaire), 1l faudrait écrire : 
dr) = À Air) utlku, r). 
t 


Chacune des F4 devrait alors satisfaire à l’équa- 
tion (13) contenant le w correspondant. 

Revenons à (12), substituons V(r) dans (7) et 
utilisons l’équation satisfaite par u(K, r), déduite 
de (6), nous obtenons : 


(13) 


ri ri] , h? 9 1 
HD Fr) :- [E REX + iä] F(r). 


Nous considérons le deuxième et le troisième 
terme du premier membre comme une perturbation 
que nous négligerons. Si les w(k, r) sont très peu 
ondulées avec r, yu — 0, cette approximation est 
permise et le choix de k,, est indifférent ; si cela a 
lieu même pour la fonction d’onde du bas de la 
bande de conductibilité, nous pouvons aussi bien 
faire 4, = 0 et utiliser (5) au lieu de (12). Nous 
avons alors aussitôt l'équation de Wannier avec la 
masse normale : 


2 
AE E Vi Ft) = EF(r). 


Cela correspond au cas où la masse effective est 
précisément égale à la masse normale [3]. 
Lorsque ceci n’a pas lieu, la perturbation est 
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minime si *, est choisi égal au moment moyen. 
Dans le cas : 


V; = constante Ir 


eik,,.r 
ces deux termes se détruisent et (13) se ramène à : 
ME Eilk). 


On a retrouvé (9). Dans le cas d’un potentiel 
quelconque ces deux termes ne se détruisent pas, 
mais la perturbation du premier ordre qu’ils intro- 
duisent sur l’énergie est : 

ñ°? Vu LÉPRE VU 

css Lt D) RUN ES 

fr 9 F dr + 5 2ükn JF LE F dr. 

Sortons le terme yu/u qui a la même valeur 
dans chaque cellule, cette perturbation est nulle si : 


[r*v Fdr — ükn JF F dr 


ce qui traduit bien que k, est le moment moyen 
(à 3 dimensions, la démonstration est plus com- 
pliquée). 

Avec ce choix de k,, et si de plus u(kh», r) est 
faiblement ondulée, l’équation : 


n° al 5 1 à ; ñ? Ë 2 
nes DUR AP [2 — Enlkn) +5 | E(r) (14) 


sera vraisemblablement une bonne approximation. 
Elle se réduira à l'équation de Wannier avec la 
masse normale : 


EP) + M Fi) ie Fr) (2) 


2m 


s’il est permis de poser : 
En (Km) = Ge ko 


c’est-à-dire si la courbe e(k) se confond au niveau 
de k — k, avec la parabole des électrons libres 
ayant pour origine le bas de la bande de conduc- 
tibihté (4 = 0, & — 0). 


En conclusion, nous avons signalé au $ I l’accord 
avec l’expérience de l’équation (2) (équation du 
type Wannier, mais renfermant la masse électro- 
nique normale et non la masse effective au bas de 
la bande de conductibilité) pour le calcul des 
niveaux localisés : 

a) suffisamment profonds au-dessous de cette 
bande, mais dont la profondeur (quelques dixièmes 
d’électronvolt) reste néanmoins encore faible devant 
la largeur de la bande interdite ; 

b) pas trop fortement localisés. | 

Si cet accord n’est pas le résultat d’une coimci- 
dence fortuite, il nous faudrait admettre : 

a) qu'à une hauteur de quelques dixièmes 
d’électronvolt dans la bande de conductibilité le 
comportement des électrons est voisin de celui des 
électrons libres, ceci étant supposé avoir lieu même 
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si la masse effective au bas de la bande diffère appré- 
ciablement de la masse d’un électron libre ; 

b) que la hauteur de l’électron dans la bande 
serait assimilable à l’énergie cinétique de celui-e1. 

Ces hypothèses ont été faites au début dans les 
modèles «naïfs » de cristaux. L’étude de modèles 
plus récents nous montre que ces hypothèses ne 
sont pas réalisées d’une manière générale, mais le 
sont dans certains cas particuliers. Mais ces modèles 
plus récents sont eux-mêmes trop simples et on 
suggère que le point de vue «naïf » aurait en réalité 
une validité assez générale, surtout dans le cas où la 
bande de conductibilité est large et séparée elle- 
même de la bande de valence par une large bande 
interdite. En effet, plus la masse effective au bas 
de la bande de conductibilité sera voisine de la 
masse électronique normale, plus ’mtervalle d’éner- 
gies où cette approximation est valable sera grand. 


Appendice : 


Quelques propriétés du réseau de Kronig et 
Penney. — 1. Considérons le réseau à une dimen- 
sion de Kronig et Penney [14] : pointes de potentiel 
en à pour æ — 0, + a, + 24, etc. On posera eV 
potentiel dans le réseau non perturbé : 

J à + co 
OR RS ENS 
ee ma | à 


j==69 


d(a — ja) 


c étant le paramètre ainsi désigné dans Seitz [14]. 
On obtient une fonction de Bloch de nombre 
d'onde À : 
Œ(k, æ) — uk, x) eike 
u(k, & + a) = u(k, x) 
correspondant à lénergie £ pour : . 


C : 
COS & & + — sin &œ « — os ka (15) 
œ 


a = V2mE/n. 


Les limites de bandes sont données par ka —n 7, 
doute à 


COS ay & + sin &y a@ — (— 1}? 


Cr 


Par développement limité de (15), on obtient. 


aisément la masse effective au bas de la nieme 
bande, définie par 


2 — 1 
MI 01 ( ) 
dk? 
m* : snæa, € Sin dy & 
D t É COS an a)|. 
m An & ia) Nate 


(16) 


On voit que »* au bas d’une bande tend vers 
zéro quand cette bande est élevée. 
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La masse effective dans la partie centrale d’une 
bande élevée tend par contre vers la masse de 
l’électron libre (*). Pour « 4 grand, on peut en effet 
négliger dans (15) le terme € sin « a/x a devant 
cos « a (sauf quand ce dernier est très petit c’est-à- 
dire au voisinage de la limite de bande), ce qui 
revient à faire # = à. 

2. Considérons le niveau «d’impureté » obtenu en 
abaissant de W le potentiel d’une ‘cellule 
(V;, = — W}). Ce cas a été traité rigoureusement 
par Seitz [14]. En recherchant la solution qui 
s’amortit de exp (— x a) dans les cellules succes- 
sives, on obtient la condition de raccordement : 


Yi Ya 


EURE Per + RE à —U4& 
Bo Sn EE COS me He 
ou 
NES 
T'eot = sin œ'@ + COS mi — + CT 
(4 


à joindre à (15), ce qui s'écrit alors : 
C : % 
COS &œ a + — sin «a = + cosh ua 
x & 
VomE 


VEntE + M) 
x — A, 
h À h 


3. Comparons la solution rigoureuse ainsi 
obtenue à celle donnée par l’équation de Wanmier. 
On sait qu’un trou de potentiel de profondeur W 


- TES 


FiG. 5. — Niveaux d'énergie d’un troù de potentiel. 

Notations : W profondeur du trou dont le sommet Æ, 
se ra au bas de la bande de conductibilité, Æ niveau 
ocalisé. 


et de largeur à (fig. 5) possède une série de niveaux 
donnés par 


ee FA 
ts TE VIm à = \/£ ou bien 4/2. (17) 


A) Lorsque le trou est très peu profond, (17) 
devient : 
DIT 


ee TU Ce mW?, 


(*) La fonction de Bloch du bas de la bande k = 0 
est 9: : COS x (2 + tandis qu’on vérifie que les fonc- 


tions de Bloch de la partie centrale des bandes élevées sont 
faiblement ondulées (4 = const.), 


L 
{ 
4 
i 
: 


noté ti dit mnt ft Mu ÉEÉ téRe de SEL à 


. 
| 
| 
| 
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Un calcul un peu long, mais élémentaire, montre 
que la solution rigoureuse du $ 2 devient bien 


alors : 


DÉTIUE 1 
E— m* W=? 


h? 


F16. 6. — Réseau de Kronig et Penney dont une cellule 
est profondément altérée. Le niveau £ est très localisé et 
son énergie s’identifie à celle d’un électron de masse 
ordinaire dans un puits de potentiel très profond. 


c étant la profondeur du niveau localisé £ au- 
dessous du bas de la bande de conductibilité £;, 
et m* la masse effective donnée par (16). 
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B) Considérons le cas où W est très grand. 


L'énergie du niveau localisé devient alors néga- 


tive — |E| (fig. 6). 
On posera B = V/2m|E| [h, et les formules du $ 2 
deviennent : 


< Fe sinh $ a + cosh B a — e—ua 
cosh 8 a + ga sinh Ba— + cosh ua. 


Si |£| elle-même est très grande, l'amortissement 
de la fonction d’onde du niveau localisé hors de la 
cellule altérée est très rapide (exp (— uw a) — 0), 
il vient sensiblement : 


pe ES er 
Lg L* (ou cot L*) co 


HAE 
Le. "2 Da SR 
8ma? 


ce qui est bien l’expression des niveaux du trou de 
potentiel de profondeur infinie, avec la masse ordi- 
naire 7n et non m*. 

On retrouve ainsi, qualitativement, certaines des 
propriétés générales admises ci-dessus. 


Manuscrit reçu le 12 mai 1958. 
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LES MODÈLES DE FORCES INTERATOMIQUES DANS LES SOLIDES 


Par Y. LE CORRE, 


Laboratoire de Physique Théorique, Collège de France, Paris. 


Résumé. — Les différents modèles de forces interatomiques dans les solides sont discutés. 
On donne l’expression générale des constantes de forces CÈ compatibles avec les théories de 


Born et de Laval. Les coefficients a ne sont pas symétriques en « et 8. L'écart aux relations de 
Cauchy n’est pas dû aux forces non centrales mais au fait que la périodicité cristalline est per- 


turbée près de la surface. 


Abstract. — The different types of interatomic forces in solids are reviewed and discussed. 
The general form of the force constants GE. according to Born’s and Laval’s theories, is given. 


The coefficients CS are not «, B symmetric. The deviations from Cauchy’s relations are caused 
not by the non-central forces, but by the fact that the crystalline periodicity is disturbed near the 


surface. 


Nous passerons en revue les différents modèles de 
forces interatomiques et examinerons les conclu- 
sions qu’apportent chacun de ces modèles à la 
théorie de l’élasticité (relations de Cauchy, théorie 
de Laval). 


1. Forces dérivant d’une énergie radiale. —- Nous 
dirons que l’énergie du solide est radiale si elle ne 
dépend que des distances relatives moyennes des 
ions entre eux : 


æ = (ro) ryn —= distance de l’ion g à l’ion . 


Un cristal dont le champ de forces dérive d’une 
énergie radiale satisfait à la théorie classique de 
l’élasticité. 


1. CAS GÉNÉRAL. — a) Les variables indépen- 
dantes. — Considérons un système de V atomes ou 
ions, soit g un ion particulier, ses coordonnées 
seront notées : 


2% (x, B, y = 1, 2, 3) 


la distance relative de déux atomes gh est notée 
Ton — Th: 

Autour d’une position d’équilibre, la configu- 
guration du système est déterminée par 3-6 
variables indépendantes r,, définies comme suit : 

On utilisera les coordonnées tripolaires en 
prenant pour pôles 3 atomes e ou f — 1, 2, 3, la 
position d’un atome g légèrement déplacé par 
rapport à la position d'équilibre est déterminée par: 
lois l'o2s l'as: 

On constate que les 3N-6 variables indépen- 
dantes sont tous les r,, ou r,, : 
pou h'= 1,2, ONU LT Or RrchES} 


b) Expression des constantes de forces atomiques. 


La différentielle première de l’énergie est nulle à 
l'équilibre : 


du ue Me 
D Sn) Le 28) (328 — Bag) = 0 (1) 
donc : 
ù w dw 
—— — (). 2 
Dry oUréd @) 


Il en résulte que la différentielle seconde s’éerit : 


C7 
Bp=2 D YD = ec 


g>ea h>1.8 d(r*e) d(r?) 


(x% — 2%) (28 _ a%) (Sa% — da$) (524 — Bab). (3) 


Posons : 
LR h, ==4l, 2; 4 
£ d?w si g> En RE 
ce MST ENT D ‘ 
+ d(rÿk) d(rÿ1) k, l == ts 2, 3 
nl gi, — : d . 
0 si R PINOT 
ou g<kou k << 
On note que : 
Txnt ES Tu, or (5) 


et d’après (1) 


9 1 a 7 , 1 
Sw => ZE (Ten + Lori + Leon + Ty) 
2 gkx hf 


ou encore : 


D? w — na Z C2 8aÿ Br$ (7) 
2 gahf 
avec 
h 
CR = — 


! * ! ' ” \ h l 
> (Lou + Toxin + Tront + Lroun) (x — 24) (28 — x). 


(8) 
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On note que, d’après lexpression (8) de C#£ e 
d’après (5) 


CR = Che (9) 
D cn = (10) 
>) CE a - -È CÀ ab. (11) 


Les relations 11 ont été établies par Born en 
considérant une rotation rigide du solide. 
En général, C%£ n’est pas symétrique en « et 8 


m £ Cêa. (12) 
c) Expression des forces de rappel. — Soit F£ la 


composante de la force exercée sur l’atome g, on 
doit avoir : 


dd — D RENE 
gx 
c’est-à-dire : 


bin 007 ah sh. 
F£ = sas Pa Cu 5 x (13) 


La constante CŸ£ a une signification simple : elle 
donne la valeur de la force à laquelle est soumis 
Patome g dans la direction « si on déplace l’atome À 
dans la direction £. 


d) Forces compatibles avec le déplacement d’un 
atome. — Déplaçons l’atome gseul, s’il est soumis 
à la force 


fe — Co D af 


cependant que les autres atomes k (h  g) sont 
soumis à la force 


fa — Caë 
On constate, d’après (10) que-: 

9 Te , 

fa + ZE fa = 0 (13) 


Ê 5 af. 


et, d’après (11) 
fa 2% — 4 2x si 2, (fa — j} 2) At) 


c’est-à-dire : 


Moment f9 + X Moment f# — 0. (15) 
h£9 


Les égalités (14) et (15) signifient que les forces 
compatibles avec le déplacement d’un seul atome 
ont une résultante nulle et un moment résultant 
nul. D’après le principe d’action et réaction ces 
forces sont égales et directement opposées à celles 
qu'il faut appliquer aux différents atomes pour 
déplacer l’atome g seul. 

Ainsi les forces qu’il faut appliquer pour déplacer 
un atome ont une résultante nulle et un moment 
résultant nul. 
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e) Constantes élastiques des cristaux déformés 
dans une affinité. — Dans un cristal, nous représen- 
terons union par deux indices "7 affectés, le premier 
à la maille où il est situé, le deuxième à sa position 
dans la maille. » est l’affixe du centre de la maille, 
c’est uné combinaison linéaire des translations 
élémentaires. » + j est l’affixe de l'ion nj lui- 
même. 

La différentielle seconde de l’énergie s’écrira : 

Mon es ŒIL Sr) à 2. (16) 


2 mi) 
20 


Une déformation par affinité est définie par les 
composantes du tenseur des déformations {8 


è 2 =. © (me + J8) tag. (17) 


L'énergie de déformation est : 


D? w — —— Es C 5k | (my + jy) (ps + ks) tay tps. 
mjay  %P 
CUMCN) 


Il vient pour un £ristal de volume V : 


OT Us tone 65 Lay !BS 
2 œy,B3 
avec 
PL ce 
M oy,85 HR he Er ps (18) 
mi 
pE 


d’après la relation (11), les coefficients d’élasticité 
sont symétriques en « et y, en 8 et à, en «y et Bñ. 
Il n'existe que 21 coefficients ON, 45 indépen- 
dants. On pouvait d’ailleurs tirer directement 
de (3) : 
4 d? w 


sie TE 
rer ur a>eu h>1,6 07%) (réf) 


(29 — 29) (ag — 29) (2ÿ — 28) (28 — 26). 
f) Constantes élastiques, cas général. — En 


général, les sous-réseaux subissent des translations 
relatives : 


m 


Ô a j — a, + 2 re 
D’après (10) la densité d’énergie peut se mettre 
sous la forme : 


— Er + > Mylay) (EG — E) 
+ > PS 88) |: 


Si la maille contient g atomes, il existe g-1 dépla- 
cements vectoriels relatifs indépendants # — El. 
Ces déplacements sont obtenus en rendant l’énergie 
de déformation minimale : 


mo . 
Z CHE — & + D pe ts) = 0. (19) 
"r 
1100 
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Appelons S la matrice à 3g-3 lignes et colonnes, 
inverse de la matrice dont les éléments sont : 
—XC% j44,k%1 


c’est-à-dire que : 


DD CIE = 1. (20) 
ja,k@ mp 
Orrar 
de (21) 
Il vient : ' 
-h My 
Pete ZE 44 CE pe 5 (22) 


mix 
pes 


et, portant dans l’expression de l’énergie : 


F : 
She F5 2 May.ps 
— «y, Rà 
1 
To Re my Ci ê è à Se Cikpsl] luyles 


Mi,phk nie, den 


Les coefficients d’élasticité sont donc : 


: ; AT: 
Moss = Morts + May.s5 (23) 


avec : 


ph à 1 + mn 4 
My: CRE >> * > nty C ee Lee en C k ps 
} Mmji,pk Nie,qen 


ou encore : 


1 + 
D 


û NE B5 — y c : 
mi,pk mie,qer 


MAP TE 6 AD 
(my — ny) ( 25 Sen Ceé (gs — ps) 
(24) 


Les relations (11) imposent à nouveau aux coef- 
ficients AU 45 la symétrie en « et y, en f et à, 
en «y et BS. Cest également la symétrie des coef- 
ficients d’élasticité 9T,,g5s. Les tensions sont 
exprimées en fonction des déformations par leurs 
composantes : 


= 2 My.ps le: 


Le tenseur des efforts est donc symétrique en « 
et y 

Les 21 coefficients OT, #5 sont les coefficients de 
Voigt. 


2. LES APPROXIMATIONS. — a) Forces dérivant 
d'énergies biatomiques. — Certains auteurs admet- 
tent implicitement que l'énergie peut se mettre 
sous la forme 


C, (5 29 2% — 5x4) (3 28 — 5 x$). (25) 
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On voit qu’il est nécessaire que les termes en 
dx das et dxf dax soient égaux : 


= As CB (26) 


la condition est es car, compte tenu des 


relations générales (9) et (10) 
8 — ; (9) 
D EIRE (10) 


On peut écrire successivement : 
y — AS 5 252 Ca B(S 28 — 5 2%) 
—;2 5 a > CRE(5 x8 = D) 


: | 
=> D CÉ(S 2% — 5 àù) (5 28 — 5 28). 


4 ghaf 


L'énergie de déformation s’exprime done comme 
la somme d’énergies relatives biatomiques 


avec : 
02 gygh = F2 jphg == = , 2 CS 2% Lx — 0 x!) (5 28 0 24) 
et 


, gh d? #w9h 
C5 af — 


7 (27) 
das — 1) (af — 24) 

On peut alors donner un sens physique à la force 
relative f*", qui n’est fonction que du déplacement 
relatif des atomes get À par rapport au cristal non 
déformé : 


ah dwuh LE ct h 
= — —————— C 2 ACET TG — ù TL k: 
ee dl 00 ES 


On note que f#" + jh" — 0. 
La force globale à laquelle est soumis l’atome g 
est : F9 =2 {g". Notons qu'il est toujours possible 


d’ effectuer. ‘cette décomposition de la force glo- 
bale F% : 


pe -È AMEN Aa —È CS 28 — 5 a). 

Mais ce n’est pas la seule décomposition possible, 
f£" apparaît comme un artifice de calcul comme le 
champ élémentaire de Laplace en électroma- 
gnétisme ; la force relative {g" n’a un sens physique 
que si elle dérive d’une énergie, c’est-à-dire si C2g 
est symétrique en « et £. 

En élasticité cristalline, IÙ,,#5 peut se mettre 
sous la forme simple : 


Mes = 
1 , mp ‘ 
Eve ae (my — py + Jy—ky) (ms — ps + 15 — 3) 


(28) 


No 8-9 


Inversement, on ne peut généralement 

mettre AN, ss sous cette forme que si Ce est 

E «x 
symétrique en «x et GB. Les coefficients AU, 5 
deviennent done symétriques par rapport à leur 
quatre indices. C’est-à-dire que si chaque atome 
est en un centre de symétrie, les relations de Cauchy 
sont satisfaites. 

b) Cas d’un cristal où le champ de forces a un 
rayon d'action sensible. — Nous dirons que le champ 
de forces a un rayon d’ation sensible si Ci décroit 

: re “ . Œ 
vite avec la distance des ions m»7 et pk, et ne 
dépend que des voisins situés également à courte 
distance de "Jet pk. 
EE mp : . 
On pourra écrire que C Le est invariant dans une 
œ 


. translation propre du cristal et ne dépend que de la 
différence M — p : 


mp MD 
FN es 


De plus ce coefficient est soumis aux éléments de 
symétrie qui conservent ou permutent les 
atomes 7717 et pk. 

Il en est toujours ainsi sauf quand le champ 
coulombien introduit certaines séries semi-conver- 
ventes. Ainsi, dans un réseau de points chargés dis- 
tribués comme dans un cristal de CINa, la somme 


M —D 
DO,» 
(4 


m—p % 


ik 


prend des valeurs différentes selon que l’on somme 
une sphère ou sur une plaque normale à Ox; 


= MP 
DS RUTENE VE) 
œ 


1] 
mp 


sphère 
> Ci à 167 
D D 77. 
. #j <k 
AD 3 U 5 
plaque 
D co D 147 
DCE Ci Ts 
mp 1 Many AN IUT S 
plaque plaque 


v est le volume du motif (maille rhomboédrique), 
z; et z, sont les charges des points 7 et k. Ce résultat 
est en accord avec les diverses expressions du 
champ dépolarisant. Cependant dans les réseaux de 
points chargés, c’est la seule difficulté que l’on 
rencontre. Les conséquence de cette semi-conver- 
gence se limitent à l’expression de la fréquence 
optique de vibration du réseau. Les séries con- 
tenues dans l'expression des constantes élastiques 
sont absolument convergentes. Il n’en est pas de 
même dans les métaux : le déplacement d’un ion 
métallique »m perturbe lion q qui réagit à son 
tour sur lion p. Cet effet n’est important que pour 
les ions situés près de la surface. 
L'expression du coefficient 


Ca8 
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dépend en grande part des ions q éloignées, son 
expression ne dépend pas que de la distance rela- 
tive m-p. Le coeflicient C#Z est différent pour deux 
couples d’ions situés à la même distance relative 
mais en des positions différentes par rapport à la 
surface. 

Si le champ de force a un rayon d’action sensible, 
de sorte que les écarts à la périodicité cristalline 
près de la surface soient sans importance, l’expres- 
sion de JN,.,.s se simplifie : 


ME : Ci k 
SRE = 7j Emyps (30) 
1 MD 
SC Ci km ep) p 
V mi.pk X (9 ù Py) . 
il M—D 
=. > ap ZE Ci L ps 
V M—D,i 
1 M —D 
= — Cam NES 
2, : à EU — pr) (ms — ps) 
1 ( M—P 2 3 À) 
= — à RCTETE — ES 
2, > Ci £ | Ci . (ny — p4) (ms — ps) 


V est le volume total du cristal, v le volume du 
motif. 

On constate que AN, 45 devient symétrique par 
rapport à ses quatre indices et vérifie les relations 
de Cauchy. Seule la partie symétrique des cons- 
tantes de rappel intervient dans son expression. 

c) Cas d’un cristal à un alome par motif où le 
champ de forces « un rayon d'action sensible. — 
Considérons un cristal contenant un seul atome 
par motif. Le point milieu de deux atomes g et À 
est centre de symétrie de tout le cristal si on ne 
tient pas compte de ses limites (cas où les forces 
entre atomes possèdent un rayon d’action sensible). 

Déplaçons l’atome À seul de la quantité 32% dans 
la direction 8, on a ; 


SUN 
É = Ca ù &p. 


Si on déplace atome g de la même quantité, 
mais dans la direction opposée, il vient : 


fa — — Cf 3 a 
mais le centre de symétrie impose 
fa = — fa 
donc : 
ch = C8. 


Cette symétrie en g et k, donc en «x et & permet 
d'écrire l’énergie sous forme d’une somme d’éner- 
gies biatomiques. 

d) Cas d’un cristal où le champ de force a un 
rayon d'action sensible et où chaque atome est en un 
centre de symétrie. — Appelons à nouveau 7] 
(ou pk) l’atome d’espèce 7 (ou k) engagé dans la 
maille » (ou p). Son affixe est »m + 7 (ou p + k). 
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Si chacun est centre de symétrie, on a : 
] —Rk =T.}2 
où T'est une translation élémentaire. 
On peut démontrer, comme précédemment, que : 
ner p—T—Mm 


C5 Ps % 


CAC) œ f 


(31) 
il en résulte que : 


S CE 
œ 


mi,Dle 


M—D 


k + C 5 }) My PS 


6 
2 


D—mMm ml 
é C5 + ch 1] My ps 
mi,pk Ÿ 


(Q 


c’est-à-dire que A, es est symétrique en « et $ 
et peut s’écrire : 


oo +0) 
oy,ps — AD 2,2 fer & FC é « 
(my = py LI; — : Yl' (ms ps + j5— Ks) (32) 


MI es étant nul, on retrouve ainsi directement 
les relations de Cauchy. 

e) Cas des cristaux métalliques cubiques. — Ces 
métaux ont un seul atome g (ou k) par motif. Nous 
écrivons : 

Cas = CE 6 + C'a 

C% 8% est le terme périodique que l’on trouve loin 
de la surface. Ce est, au contraire, un terme qui 
n'existe que près de la surface et s’annule rapi- 
dement dès que l’on a pénétré de quelques dis- 
tances réticulaires dans le cristal. 

Ce terme apporte à l'énergie une contribution 
créée par le travail des forces subies par les atomes 
superficiels : 


(33) 


Soient » la direction de la normale à la surface 
aux environs des atomes g et h et t, et {, deux 
directions parallèles à cette surface. La symétrie 
locale du tenseur C'%$ impose : 


vuh 9h : 
C Ni — C'nts = 0 ; 


vgh v Ch 


CE C 2 Can Ci = Cat: 
Il vient : 
à? a” 
—52 C'An8 a4 3 an + C'Pe (8 af, + Sat) (S 2h + Sa). 


Le deuxième terme apporte une contribution 
globale nulle car on peut toujours trouver, de part 
et d’autre du centre de l’élément de surface consi- 
dérée, des déplacements ôx? et Daÿ égaux et opposés. 

Dans le premier terme dx? et 52} sont tous deux 
sensiblement égaux au déplacement normal de la 
surface dx et il reste : 


1 
D? @" — —— 
2 
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Ce ne dépend pas de l’orientation de la surface 
puisque le cristal est cubique. 


> C'?} est proportionnel au carré du nombre des 


atomes situés près de la surface. Il vient donc, en 
admettant que 9? w’ est du même ordre de gran- 
deur que 3? w : 


LEA 
TS Œ Ô LA AS)" 
ou encore : 


LINK 
: op Vie. 

La partie non périodique du champ de forces 
apporte donc une contribution proportionnelle au 
carré de la variation de volume 3. 

SV = Vitis + too + ts). 


Les constantes élastiques seront donc : 


Doy.p8 = : Mer gs + Mor.0s (35) 
avec : 
Moy.B = >= Ra CE der = hy) (g8— hs) (36) 
May88 — _ Day Vp5- (37) 
On retrouve les relations de Fuchs : 
Tite O0 REMOTE (38) 
DÉC TRES (39) 
La relation de Cauchy : 
SRE Fe UE 
permet de connaître expérimentalement X : 
K — 90 (A, — Ni). (40) 
f) Les forces centrales de Born. — Un cas parti- 


He de décomposition en énergies biatomiques est 
celui où chaque 3?w°* ne dépend que de rm. Le 
champ de forces qui en résulte est celui des forces 
centrales de Born. 


— =; » 9h (r0p) (41) 
Il vient : 

d 1 , dw0h h h 
dw — J 2 2 EE Ge a ÉFEES a*) (42) 
et, puisque le cristal est en équilibre, 

dwoh n DNS 

ñ d(r5n) Gé — 2e) = 0 (43) 
et, 
1 d? LEE 

a SCT SP uS 
tu = À RE [AE 20 (af— cp DR + 2 DE] 


(5 2% — 5 2ù) (5x8 — 5 r8) 
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(34) : 


| 
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ou encore : 


V9, 1 NE Y) NS 
D — —> À Cnafo as D 28 (44) 
2 go, hf 
avec : 
, gh { h 9 n\ d? 0h dsl 
Cpog — 428 — ta) (xf — 28 ss 2 45 
; à nr EN 6) d(rün) | Le d(rËn) +) 
et, par convention : 
; , gh 
De 2er (46) 


| On vérifierait aisément sur ce cas particulier que 
les relations de Cauchy sont vérifiées si chaque 
atome est un centre de symétrie car AN}, prend 
la forme : 


y.e5 (47) 
fl D? 9h 
— s Ps AGE TE a}) (x8 — 28) (aÿ — x4) (x$ — a%) ST 
g) Les énergies triatomiques radiales. — Il est 


toujours possible d’écrire l’énergie de déformation 
sous forme d’une somme triatomiques 


dun > D? gefi 


“ef 


(48) 


où 3° w‘' ne dépend que des déplacements relatifs 
des trois ions e, f et t repérés par rapport au cristal 
libre de toute contrainte. 
Un cas particulier est celui où 9? °° ne dépend 
. que des distances relatives r,;, r;; et rs. 


D? wweli — D? welifré, rh, rie) 


Il vient : 
11 h\ (.E l dei 
D rairnste re 70 AE æ 
St aai = 2 [5 &(xg — za). (x6 — x8) NERE 
CARAETA A 
Da Lo 
+25 _. Bag Da. (49) 


Cette approximation représente une amélio- 
ration sur les forces centrales de Born mais, ainsi 
que l’étude générale l’a montré, ne permet pas 
d'expliquer les écarts aux lois de Cauchy, dans les 
cristaux où le champ de force a un rayon d’action 
sensible. 

On pourrait généraliser en considérant des 
énergies p-atomiques qui ne soient fonction que des 
distances relatives des p atomes entre eux. 


II. Forces ne dérivant pas d’une énergie radiale. 
— Certains auteurs ont utilisé des champs de force 
qui n’obéissent pas aux relations de Born : 

ZE Ca ay — Z Carat. (1) 

L'énergie du solide n’est plus fonction des seules 

distances relatives des atomes. Si on applique des 
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contraintes au solide, la déformation qui en résulte 
doit être rapportée au même référentiel que les 
contraintes appliquées. Ce référentiel est déter- 
miné avant toute déformation dans le solide libre 
de toute contrainte. 


1. CAS GÉNÉRAL. — Un ion dans un solide n’est 
pas un point sphérique. On le représente souvent 
par un corps central d’où partent des bras de 
liaison dirigés vers ses proches voisins. Cette confi- 
guration lui confère un moment multipolaire 
(moment dipolaire dans le cas des cristaux ferro- 
ou antiferro-électriques). Une énergie est nécessaire 
pour faire subir une rotation «, à ce moment 
multipolaire dans le champ créé par les proches 
Voisins. 

L'énergie interne du solide est donc fonction non 
seulement des distances relatives entre ions, mais 
encore des orientations relatives des ions entre 


eUX : 
D Ton 08) 


A l'équilibre les orientations sont telles que : 


d w 
4 

Dans toute déformation les orientations sont 
déterminées par la condition de minima : 


DO HE 
USE 


La déformation doit être repérée dans le même 
système d’axes que les tensions appliquées. La 
différentielle seconde de l’énergie se met sous la 
forme : 

52 pe 1 » 9h & 9 R 
W = ——> Cap d to d xp. 


2 ga,h6 


(50) 


Le coefficient C% est symétrique en gx et AB car: 


CR Ris (51) 
dr dr 
En particulier : 
CE CE (52) 
mais, en général : 
Cas # Cha. (53) 


La forme exercée sur l’atome g est : 


Il importe que cette force soit nulle si le solide 


subit une translation d'ensemble. Il en résulte : 
> C0. (55) 
h 


Cette relation impose l’invariance de l’énergie 
dans une translation d'ensemble du solide. Elle est 


710 


valable en statique si on n’exerce aucun effort pour 
retenir le plasma électronique, mais ne tient pas en 
dynamique où la force d'inertie de ce plasma entre 
en compte. 

Déplaçons l'atome g seul, il est soumis à la force 


fa — Caë 2 
cependant que les autres atomes X sont SOUMIS AUX 


forces 
fa — Cë af 


on constate que : 


Le 2,f = 0 à (56) 


mais le couple résultant n’est pas nul, il est équi- 
libré par le couple en volume qui crée un moment : 


(2x9 Pre À, RME Pa) 


K, PP 
0) 


(52) 


LE be 


(permutation «, 6, + paire). 
Ainsi les forces {et Y 2 sont égales et opposées 
RF 


g 
“mais ne sont pas directement opposées. L'énergie est 
inchangée dans une rotation d’ensemble car le 
travail du couple + annule celui des forces f*. 
L’élasticité cristalline se traite comme dans Je cas 
des énergies radiales ; 


Æ 1 PHSPTT A 
M y,65 = a — NT ,.65 (58) 
<T mp : + 
My.ps = EC 7h (my + jy) (ps + ks) 
. mpik 
1 mp É 
——+ 2 Cikmyps (59) 


Mass = F2 (my + jy) CÿE S2 CIE (ps + Ks). (60) 

Les coefficients IX, et As ne présentent 
pas, en général, d’auire symétne que celle en &y 

83. En particulier, la relation (11) n’étant plus 
valable, OT, +5 est différent de ML, ;s et le tenseur 
des efforts 1°,, n’est pas symétrique en & et y. Le 
nombre des constantes élastiques indépendantes 
s’élève à 45 pour un cristal triclinique. Si chaque 
atome est en un centre de symétrie, AMLes est 
encore nul. I] nous faut reprendre l’étude des urvers 
cas particulier sans prendre en compte les éga- 
lités (8) et (11). 


2. CAS PARTICULIER. — Forces dérivant d'énergies 
biatomiques. — Le raisonnement effectué dans le 
cas des énergies radiales est toujours valable. 
L'énergie de déformation peut se ramener à une 
somme d’énergies biatomiques si 


CB — Cha 


Il vient 


2 w —— ER 2? y9h 
Zah 
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N° 8-9 : 
avec 


= 5 r%)} (5 28 = Safe (61) 

On a alors le droit de définir la force relative f? 
exercée sur l’atome g par le déplacement relatai 
(3x? — 5x") et qui ne dépend que de ce déplacement 
relatif : 


per PTE (62) 
Ta 


ZX CS 28 — 5) 
5 36 z Sa ss CD 


On note que - 


et 


(63) 


+ je = 0. e 


La symétrie en x et B abaisse le nombre des 
coefficients JT, e5 à 36 et permet de les écrire 


May.e8 


1 mp ; re 
= — PE à © + —— _— 4: Fe 5 
2V 5,0 CS (ny + pr — ki) (ms + Îs — ps — ka). 


(64) 


Si chaque atome est centre de symétrie, on doit . 
trouver les 9 relations de Laval entre les coefficients 
élastiques à deux indices : 


ER PE PP 
Is A y JL; a Is Den su TT 46 (65) 
Moy = Moss Mag — As y — gs: 


Ces relations remplacent les relations de Cauchy - 
de la théorie classique qui suppose, entre autres 
égalités, Jo —= A pe, 


Champ de forces ayant un rayon d'action sensible. 
— Le raisonnement déjà donné permet d’écrire : 


5) RS PERTE HS 


EU m—p.5k 


: EC b 4 (m1, — P+) (ms => ps). 
(66) 


Ces coefficients sont symétriques en x et fB ainsi 
qu’en y et à ; ils satisfont donc aux relations de 
Laval et non plus à celles de Cauchy. Si le cristal 
possède un seul atome par motif, il en résulte la 
symétrie en get h des constantes de rappel : 


Roi. FA (0 ME: 
done la symétrie en x et £. L'énergie est donc 
décomposable en énergies biatomiques. Si chaque 
atome du cristal est un centre de symétrie, on aura 
à nouveau : 


1 M—D N—P 
ait — ME = — ( 74 k+ Ci k}- 
ay, Coy,p5 4v 2. C 3 ÿ ê k) 


& 
(my + jy — py — ky)(ms + js — ps — ks). (62) 


Ces coefficients satisfont bien aux relations de 
Laval. 
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Cristaux où la structure impose les relations de 
Born loin de la surface. — 11 existe certains cas 
particuliers où les relations de Born (11) sont 
imposées loin de la surface par la structure cris- 
talline. 


mp LS VERTE 
2 CIE (Pr + k;) = 2 CiX (pe + ke). (11) 


ou encore : 
Mm—Dp m—D 

7 rs C5 E (my — p;) se 2- C5 E (reg — pe)- 

Ainsi si chaque atome est centre de symétrie 
chacune des sommes est nulle car elle doit à la 
fois changer de signe et être invariante dans une 
inversion des axes. De même, si chaque atome est 
dans un champ de symétrie tétraédrique (blende, 
diamant), les axes ternaires étant pris suivant 
(+ 1. LE 1, +1) toutes les sommes sont nulles sauf 
six, Bet y sont tous différents. 

Les symétries ternaires imposent alors : 


et les symétries quaternaires inverses : 
mp D, : 
PUR Em ps) 02 Ci E(m, — pi). 
2x RS EE 


Mais il ne suffit pas que les relations de Born 
soient vérifiées loin de la surface pour que les 
coefficients élastiques 9T, +3 deviennent symé- 
triques en 8 et à. Il en sera au contraire ainsi si les 
relations de Born étaient vérifiées pour tous les 
atomes mn, y compris ceux de la surface du cristal. 


Décomposition en énergies triatomiques. — Nous 
allons montrer qu'il est possible de décomposer 
l'énergie de déformation 3? «æ en énergies triato- 
miques 

no 52 32 ef (68) 
où 3° w‘" ne dépend que des positions relatives des 
atomes e, f et 1. ? 

On peut écrire : 


1 h Ç 9 
DD CR OS DE 
2 gxhe LA DR 5 
1 CR +. CR g as NIUE SE: 
= Rs (5x2 5 22)(5 26 — 5x8) 
4g9hax@ 2 7 
AC CR 
LZ- xp Ex y 9 5 _R 
> 2, 5 Ô T4 Ô 8. 
Posons : 
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8? w°* ne dépend que des positions relatives des 
atomes get h. 

On constate que 3?w peut se mettre sous la 
forme : 
La 


ni 
GC = 
x 2Æ. (70) 


. 2 
CLASS = 


On peut écrire le deuxième terme sous la forme : 


1 : 
— D D'CÉE(S: 322) (8278 — 51%) (71) 
SRE 21 _ 2 : 
LE 


avec : 


ñ Xe Eh 103 : 3 k E = 
Dé Che— Che Cas (CES — CE). (72 


Mais il existe une infinité d’autres solutions. 


Conelusion. — Cette revue systématique des 
différents modèles de forces interatomiques nous à 
permis de préciser certains points : 

19 La théorie classique de l’élasticité suppose 
que l’énergie cristalline ne dépend que des distances 
relatives des ions entre eux (énergie radiale). Cette 
hypothèse de départ étant posée, nous avons pu 
donner l’expression générale des constantes de 
forces interatomiques C2? compatibles avec de 
modèle. Ces constantes vérifient les relations de 
Born : 


oh À gh h 
ZX CE 7 = D CB 78. 
k h 


29 La théorie de Laval réfute les relations de 
Born. Nous avons montré que la théorie de Laval 
a l'avantage de tenir compte du travail accompli 
par les couples appliqués aux atomes. La possibilité 
de couples en volume autorise la dissymétrie du 
tenseur des efforts et élève le nombre des constantes 
élastiques à 45. 

30 Les constantes de forces C24 ne sont pas, en 
général, symétriques en x et 8, même si elles déri- 
vent d’une énergie radiale. 

4° Si chaque atome est centre de symétrie et si 
le champ de forces a un rayon d'action sensible, la 
théorie de Born impose les relations de Cauchy du 
type : I,, — Né, alors que la théorie de Laval 
entraîne l'égalité : Ms — Iso. 

50 Les écarts aux relations de Cauchy dans les 
cristaux où chaque atome est en un centre de 
symétrie ne sont pas dus aux forces non-centrales 
car les forces non-centrales étudiées par Born et 
ses élèves dérivent d'énergies radiales. Les écarts 
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aux relations de Cauchy, de même que les écarts 
aux relations de Laval, sont causés en statique par 
le fait que la périodicité cristalline est perturbée 
près de la surface et, en dynamique par les forces 
d’écran créées par le plasma électronique. 

60 Si les constantes de forces C%% sont symé- 
triques en «et f, l'énergie de déformation peut être 
considérée comme la somme d’énergies biatomiques. 
Si, au contraire, elles ne sont pas symétriques en « 
et B, il est possible, d’une infinité de façons, de 
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décomposer l’énergie de déformation en énergies 
triatomiques. 

70 Le modèle le mieux adapté à l’étude des 
cristaux non métalliques consiste à considérer le 
champ de force comme la somme d’un champ cou- 
lombien et d’un champ tenñsoriel. Le champ cou- 
lombien- doit tenir compte des interactions 1ons- 
ions, ions-dipôle, dipôles-dipôles Le champ ten- 
soriel sera développé en série suivant les proches 
voisins de divers ordres. 


Manuscrit reçu le 21 mars 1958. 
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ÉQUATIONS MACROSCOPIQUES DE LA RÉSONANCE QUADRUPOLAIRE 


Par François LURÇAT, 


Laboratoire d’'Électronique et de Radioélectricité de la Faculté des Sciences de Paris (Fontenay-aux-Roses). 


Résumé. — On donne des équations de la résonance quadrupolaire analogues aux équations 
de Bloch pour la résonance magnétique, et qui permettent d'étudier les phénomènes d’induction 
et d’écho en résonance quadrupolaire. On discute le degré d’approximation des équations de Bloom, 


Hahn et Herzog. 


Abstract. — Equations are given for quadrupole resonance, analogous to Bloch’s equations for 


magnetic resonance. 
nance to be studied. 
is discussed. 


I. Introduction. — Pour expliquer la plupart 
des phénomènes de résonance quadrupolaire, il 
suffit de considérations énergétiques. L’interaction 
entre le moment quadrupolaire du noyau et le 
gradient de champ électrique (d’origine moléculaire 
ou cristalline) fait apparaître des niveaux d’énergie. 
Les phénomènes d’absorption de l’énergie d’un 
champ radiofréquence et de relaxation énergétique 
s’expliquent par les transitions des noyaux entre 
ces niveaux. 

On connaît, cependant, un certain nombre de 
phénomènes de résonance quadrupolaire que les 
considérations énergétiques ne suffisent pas à 
expliquer. Il s’agit essentiellement de l’induc- 
tion [1], et en particulier des échos [2], et de l’effet 
Bloch-Siegert [3]. Ces aspects de la résonance qua- 
drupolaire, comme leurs homologues en résonance 
magnétique, font intervenir des propriétés qui ne 
sont plus «scalaires », mais où la cohérence entre les 
spins joue au contraire un rôle essentiel. En réso- 
nance quadrupolaire, on peut obtenir des équations 
analogues à celles de Bloch, bien que plus com- 
pliquées [4], [5]. 

Nous établirons et discuterons ces équations dans 
les paragraphes 2 à 5 ; nous résumerons les résul- 
tats obtenus, ainsi que leurs principales consé- 
quences physiques, dans le paragraphe 6. 


2. Méthode d’écriture et d’intégration des équa- 
ions. — Les équations de Bloch peuvent 
s’obtenir [6] à l’aide de l’équation générale : 


4a< À > 


7 F 2 
À ns (de, À] > (1) 


où 4€ est l’hamiltonien d’un noyau placé dans le 
champ magnétique, et < À > la valeur moyenne 
de la grandeur À. (On néglige ici les phénomènes de 
relaxation.) Procédons de façon analogue pour la 
résonance quadrupolaire. L’hamiltonien d’un 
noyau placé dans un gradient de champ électrique 


These equations allow induction and echo phenomena in quadrupole reso- 
The degree of approximation of the equations of Bloom, Hahn and Herzog 


de révolution autour de l’axe Oz, et dans un champ 


— 
magnétique H, est 


: (2) 
. 3eqQ JR Ro 
RQ Dre — 5 14 © Me à y: 
TES ESS AE EAN IE l 
d° © Fe PE à 
Re est l'intensité du gradient de champ. 
Si on écrit l'équation (1) pour les composantes 


> 

du vecteur < I >, proportionnel à l’aimantation 
macroscopique, on constate que, contrairement au 
cas de la résonance magnétique, on obtient au 
second membre, non seulement des composantes de 
ce vecteur, mais encore des composantes du ten- 
seur (macroscopique) densité de moment quadru- 
polaire. Si on cherche alors les dérivées des compo- 
santes de ce tenseur, on obtient en général 
(si Z > T) au second membre des composantes du 
tenseur densité de moment octupolaire, et ainsi de 
suite. Tous les moments multipolaires non nuls du 
noyau, eb les tenseurs macroscopiques corres- 
pondants, interviennent donc dans les équations 
macroscopiques de la résonance quadrupolare. 
Nous utiliserons pour écrire ces équations les opé- 
rateurs tensoriels Ÿ,,, définis par les relations de 
commutation 


(IE, Yeu] = VR(k + 1) — ul +1) Yeux; (3) 
[Z>, Yrvl = Yru 


et par la condition de normalisation : 


Ye = (— 1) (1F)E, (4) 
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Les Y,, sont des combinaisons linéaires des 
composantes del : 

Yo = V2 Yi = + IE (5) 
donc leurs valeurs moyennes sont des fonctions des 
composantes de l’aimantation M : | 

NyYA< Ypo > =V2M, (6) 
NYR< Yu > = F(M;+iM) 


(N, nombre de noyaux de l’espèce considérée par 
unité de volume). 

Les Y., sont des combinaisons linéaires des 
composantes de l'opérateur moment quadru- 
polaire Dr (?) : 


Se V:Pr 0 LE 1)] 4/52. 


a(l, IE + IET,) = + (Dys + 1Dy) (7) 


XV 


1 : 
a(1+)? = 3 (Dax — Dyy) + Dyy 


(: 4 TT “ 


donc leurs valeurs moyennes sont des fonctions des 
composantes du tenseur @:x, densité de moment 
quadrupolaire : 


Ne Vu =/5 0 
INIST MO EN O) (8) 


(Qu Fe, Qu) = 1Qxy. 


Ne vEToere = 


Les Y:,, etc.., ont des significations analogues. 

Pour écrire les équations de la résonance qua- 
drupolaire, il faut calculer le commutateur de Y,, 
avec l’hamiltonien (2), donc avec Z,?: 


[T3, Yru] A T,[Z, Yyul LU Yyul Ti u(Z; Yru 4e Yru TPE 
(9) 


Le second membre de (9) se calcule à l’aide des 
propriétés des opérateurs tensoriels [8] : 


< m|Yyulm > 


= (it < Ill > PULIR ; —m, m°, u) + (10) 


V est le symbole de Racah, lié au coefficient de 


(2) Cet opérateur est l’analogue du tenseur classique 
Dix = ZE e(3x; 2e — x r?). 

L’énergie quadrupolaire s’exprime, dans le cas classique 

comme dans le cas quantique, par : 
Æ sd ® p. 
600% 

La définition du moment quadrupolaire utilisée ici est 
(à un facteur 3 près) celle d’Ayant [11]. Le tenseur macro- 
scopique Q:x, densité de moment quadrupolaire, est défini 
(avec des conventions un peu différentes) par Rosenfeld [7]. 
Ici on a utilisé un tenseur de trace nulle, qui est donc nul 
pour une distribution de charge à symétrie sphérique. 
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Clebsch-Gordan < -—m, m'|k, — u 
tion : 


>_par la rela- 


L— m, mlk, —u > 


= (—4)Eu VO LA VUIRE mn 0) 


Nous utiliserons la relation de récurrence entre 
les coefficients de Clebsch-Gordan d'indices k, & —1 
et k + 1, qui s'écrit dans le cas présent : 

— (2m —u) < — m, mlk, —4u > 

NE R +1 +u) Qi +k+2)(21—À 

= (24 + 3) (24 +1) 

<—m,m|k+1,—=p> 

TIRE + u) C1 +1 + à @I +14 

(2k + 1) (2k —1) 
<— m,mlk—1,—u>: (12) 

Les relations (10), (11), (12) permettent de cal- 
culer l'élément de matrice du second membre 
de (9), à condition que l’on connaisse l’élément de 
matrice réduit < 1 ||Y4|| Z > qui figure dans (10). 
Ce dernier se calcule à l’aide de (4) : 


Pal > NE MINIER 
AIT 


= EAN < 114 


LÉ d _ | HEIN 
= (—1)* il VI21—1+14) =, (— DUREE (13) 
mais d’autre part d’après (10): 
Ia 
= (= 4)8L < IN YEN > VUIRS = LL 


En rapprochant (13) et (14) et en calculant le 
symbole V en fonction de ses arguments [8], on 


obtient : 
(is QUE k +1)! 
a V (@1—Rk)1(2k)! 


< 1 YA > 
Finalement, en combinant les équations (10), 
(11), (12) et (15) on obtient : 


Je Yu Te Yru 15 


(15) 


(21 + k+1)(21—k+1) /kUk + uw Kw) 
2k +1 V 2(2k — 1) Paru 
UE ir. re TS ei + TT b) Yrk+iu- (16) 


L’équation d’ RE de la quantité < Y4, > 
s’écrit donc, en omettant pour simplifier les signes 
Ter 


deu L BL KE + 1) (BI —R + 1) 
GER AE 2k + 1 
V + u) (k— uw) ; 
2(2k —1) D: 
s 2(k + 1 +u) (k +1 —u) 
DRE ETC EURE 


4 NS 
me LR AA LA + u) (k—u) HT Tru tr 
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4 ceare 


RURAL) AT 


—1y AH; u Ye. (17) 

Pour des noyaux de spin Z, les opérateurs Y4, ne 
sont pas identiquement nuls pour # < 21, [u] < k. 
Il faut donc en principe écriré 47(1 + 1) équations. 
En fait, on verra plus loin que la plupart d’entre 
elles n’interviennent pas dans les résultats finaux. 

Nous intégrerons les équations en utilisant une 
méthode analogue à la méthode du champ tour- 
nant en résonance magnétique [10]. Celle-ci peut 
se résumer ainsi : on étudie d’abord le cas où le 
champ alternatif a exactement la fréquence de 
Larmor. On applique alors aux équations de Bloch 
la transformation de Larmor, qui donne des équa- 
Lions de forme analogue mais où les termes dûs au 
champ permanent son éliminés. Cette transfor- 
mation présente de plus l’avantage de rendre 
constants les termes dûs au champ alternatif, ou 
du moins ceux («termes principaux ») qui sont dûs 
à la composante circulaire tournant dans le sens de 
la précession de Larmor. Quant aux termes dûs à la 
composante circulaire tournant en sens inverse, on 
peut les négliger, leurs effets étant très petits (sauf 
dans des expériences ad hoc). Dans le cas où on 
n’est pas exactement à la résonance, on applique 
une transformation de même forme que la précé- 
dente, choisie de façon à rendre encore constants 
les termes principaux dûs au champ alternatif. Les 
termes dûs au champ permanent ne sont donc plus 
entièrement éliminés. 

Pour procéder de façon analogue en résonance 
quadrupolaire, nous écrirons d’abord les équations 
dans le cas où il n’y a pas de champ magnétique 
alternatif, et nous chercherons la transformation 
qui élimine les termes dûs au gradient de champ 


électrique. Nous appliquerons ensuite, dans le cas 


où il yaun champ alternatif, une transformation 
de même forme, mais où la fréquence de résonance 
est remplacée par celle du champ alternatif. 


3. Résonance quadrupolaire de noyaux de spin 


1 —1. — Les équations sans champ radiofré- 
quence s’écrivent : 
d 2 T7 
d Vatr = Æ 2 Q Yoti 
(18) 
a V# mit TO 
SFAUTE SRE SERSEURS ESC 


Pour trouver la transformation cherchée, on 
écrit (18) sous la forme 


d 
di 


d 
di 


(RER 7) 07 O0 (Yes + Yo) 


(19) 
(Yitr — 


Ve F 1 Qi — Yo). 
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La transformation est donc : 
Pod (Pos de Yo) e+iQot 
Pie Pas = (Pr Var) e Fi, 
Dans le cas où il y a un champ alternatif 
Hy — H; cos ot HE 0 (20) 
la transformation sera donc : 
Peel En (Par de Fa eti@rt (21) 
neue (ES te Fa) Giee 
c’est-à-dire 
ous Fe COS Oyt + A sin ot 
Por = ets sin ot + VE COS out. (22) 
(Etre You pouru Z+fi) 
Les équations (17) pour 1 —1 (et pour 
HR" =0)s'écrivent : 
dy : t = 
se NE ON — He V2Y, 
AY x t + 
de IE 10 Yi TH (V6 Yo + DV) 
dy 1 = 
ae en MeV (Van + Par (23) 
CD t — 3 
= He VE (Ya Mg ROME | 
dYo,40 es vY y 
ER ie Ha2Ÿt 


Écrivons, à l’aide de (6) et (8), ces équations en 
fonction de grandeurs réelles. Nous obtenons : 


aM, ï 
de r:08 TE Que | 
dy 
dt E Qu) } (1) 
ñ 
= 00 Me + Y Hy (= (Qzz — On) 

dy Ê 
RCE (Qu — On) = — kr He Que 
dQuz (I) 

dt (24) 
am h 

‘& 0 One + yHe M \ 
dM, 

rs = —YH; M, 

| (III) 
( ñ ñ 
D nie) 


vx | 
se 10) EVH, (= On). | 
4 / 0 


716 JOURNAL DE 

(On a écrit le système dans un ordre qui met en 
évidence le fait qu’il se décompose en trois sous- 
systèmes indépendants. On a remplacé Q,, gran- 
deur définie de façon générale par (2), par «,, fré- 
quence de résonance, égale à Q, dans le cas présent 
(1 =1)) | 

La transformation (22) s’écrit, en fonction des 
mêmes grandeurs : 


PHYSIQUE No 8-9 
1 
My —= Mé cos y l — _. où sin &;,/ 
°4 
ee h VX 
D Que — MS sin © D + (res Quz COS © 
(25) 


% LOEB 
M, = MY cos ot + TO sin ot 
(e4 


Y À 


œ 


LA. ho %x 
= — My Sin w,t + Te Où COS of: 


En portant cette transformation dans (24), nous obtenons : 


daMx ñ H;y à 

de rie &1) — Que + + 2 Le (QE — Qu) sin 20, t | 

d ñ # ñ 
HET Q%) = (60 — 01) ME HAT (QE — OH) + cos 20 1) (n 
d À ñ 
d ee (QË = G) = — 2% H, que sin 20; + . sell +. COS 201 2] 
dQr. 
SE 2 0 (1) 

: (26) 

am sh H 1 

Le —=(@ &) Ps 2 13 [ui (1 + cos 20,1) + è T° où sin 26, ] 
dMi YA +. 

= Fe 1 [ui (1 + cos 20, t) + TE Qussin 20: ( fe 
d ñ s H à ñ 
_ ee @) = Lite bi] Mr ea [ui sin 2o,t— 7 Q%{1 + cos 20,0) 

œ 

d ñ H 
di (Lo) EE T L [ur y Sin 2w,t LA ia 2 QE COS 2:07 0. : 


L’approximation analogue à celle qui consiste, en 
résonance magnétique, à négliger la composante 
circulaire tournant en sens inverse de la précession 
de Larmor, consiste ici à négliger dans (26) les 


daM* 
dt 


LEE” 
+ (0 — ©) _ Quz — 0 


termes en sin 26, { et cos 2w, t. (On peut d’ailleurs 
montrer rigoureusement que l'effet de ces termes 
est très petit [13]). On obtient ainsi le système 
suivant : 


d ñ HÉPRO EE O 
. (> Qÿe) — (&— 01) MS at 1e Quu ; Oz =D (1) 
d (rAQu— Es). vHivt 0x 
ae & 2 vz = Ù ] 
AQ%r 
dise (11) 
27 
dMy a SH F8 
Re (© ©) = où Æ LS = 0 \ 
dM; Y-H;, * 
da Me 
(III) 
d /y# à { 
sr (= Q) + (© — o:) Mÿ RS a (Et mn) =" 
d VAS = YH; 
a (ee) (Po) = 0. AE 
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On voit que l’on passe de (24) à (27) comme on 
passe, en résonance magnétique, des équations de 
Bloch dans le système fixe à celles dans le système 
tournant : on remplace simplement ©, par ©, — w:, 
et H, (=, cos w,t) par A,/2. Cependant la 
forme de la transformation (25) montre que l’expli- 
cation de ce fait n’est pas la même pour la réso- 
nance quadrupolaire que pour la résonance magné- 
tique : (25) ne représente pas un passage à un 
« système tournant », et les termes en 7/7, négligés 
ne correspondent pas à une des composantes circu- 
laires du champ alternatif rectiligne. 

Les équations du sous-système (27-ITT) ne jouent 
ici aucun rôle. Supposons en effet qu’à l'instant 
initial le système de noyaux est en équilibre ther- 
mique, soumis seulement au gradient de champ 
électrique. Les conditions initiales sont alors : 


Qax(0) 
Qix(0) = 0 


d’autre part 


= Quy(0) = — 3 Qul0) ; 


GÆk; M:(0) = 0 (28) 


Mi (0) = 


Qi(0) — Q(0), 
Par conséquent, toutes les grandeurs que con- 
tient le sous-système (III) sont nulles à l’instant 
initial, donc restent nulles par la suite. Il reste 
donc seulement trois équations utiles (27-I) ; on 
voit qu’elles sont analogues aux équations (23) de 
Bloom, Hahn et Herzog [2], ie s’écrivent avec nos 
notations : 


M;(0). 
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dM, vH, 
dé 2 


DEN 0 C= (+ m(il—m+1)—=2 


dans le cas présent). 


(Nous écrivons A, pour désigner la quantité qui 
dans les équations de Bloom, Hahn et Herzog 
s’appelle JZ,, mais n’est pas en réalité une compo- 
sante de l’aimantation.) Si on pose comme Bloom, 
Hahn et Herzog M, — 2(u cos &w, t — v sin w,t) 
on est amené à poser dans (30) : 


ie M, = 


1y# 6 nu Que 
œ 4 


2 


On obtient bien ainsi les équations (27-I), à ceci 
près qu’il faut pour cela prendre dans la deuxième 
équation (30), non pas € — 2, mais € — 4. Cette 
différence n’est pas étonnante, puisque les équa- 
tions de Bloom, Hahn et Herzog sont obtenues en 
considérant séparément le mouvement des noyaux 
de nombre quantique m > 0 et m < 0, donc ne 
peuvent être valable que pour des noyaux de spin 
demi-entier. 

Dans le système fixe, les équation de la résonance 
quadrupolaire (7 = 1) s’écrivent : 


ñ 
LE Qui = 0 


(24-I-bis) 


d 


ñ 
= (Qu) — 00 Me | &Y H; Re = 


dé 


d É h Quu _ cz) k Ha (© * Qu) = 0 de 
œ 7 


dé 


La forme de ces équations est analogue à celle 
des équations de Bloch. 


os — ©) 0 —.0 
dt 4. Résonance quadrupolaire de noyaux de spin 
F Hope Le 1M, = 0 (30) ae Fos ne 
Fes ro, VAE Yi 
d Ya #0: (? AAA _ Y, a) AVE Yade) 
ie RRENEUTe AVE CT 
d _- Aron VE Re — "TH (2 V3 Y30 + V10 Ys,42) 
‘ Es = LOVE Fois —""Ho(VAO Yoitr + V6 Yauts) (81) 
de " REA ED PH Ce) 
Le #2 HV OT P ) 
a ca — He 2V8 (Ya SP CNE 
ue ETAT Fe 
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La transformation analogue à (22), c’est-à-dire adaptée au champ 


Hy; = H; cos ot, HR 0 (32) 
est : 
x 1 ROME Le on 
DORE Yi 5 (8 cos ot + 2) + iYo41 3 Sin out + Yade (— 1 + cos w, t) 
FETE * + 7 AN EL 
M'ÉRETREREE 1Y1,41 7 Sin Out + Yo, COS ol + He ot 
é \ 
# 6V3 ; TE TS * 4 
DR an Me (— 1 + cos wyt) + Yo ts V= sin out + Ÿst: 5 2 cos ot + 3) (33) 
6 CR 
Ts Tai cos out + DR VE sin ot 
3 
Vous = ENT V5 sin œ;t—+ Yi COS 11 
Vire You DOUT NL SC UCETTO MEL: 


Nous écrirons les équations transformées sous 
forme réelle. Pour cela nous n’utiliserons pas des 
formules analogues à (6) et (8), qui feraient inter- 
venir les composantes cartésiennes du tenseur 
macroscopique « polarisation octupolaire ». Nous 
poserons simplement : 


Yu (— 1) re 1m 
EU ET 


UE 21%y (34) 
EG 21Kyy 


et des définitions analogues pour les quantités 
transformées Y#.. 

On obtient ainsi le système des équations trans- 
formées, qui se décompose en deux sous-systèmes 
dont le second ne joue aucun rôle en raison des 
conditions initiales (cf. $ 3). Nous poserons : 


3 * 1 D. 
u = 3 N Y Ce HAUT) 


1 
U — 2N F4 nK7 


d 


Bi EN yNYE 
FES VERRE 
= Nyalit 1) (35) 
1 * 
LES Y À lo9 
Fr sNY Jo 
1 * 
y a VEN Yi Es 


Le premier sous-système s'écrit avec ces nota- 
Lions : 


l , 
_ = — (© — ©.) 0 + 3y H;(— wsin ©,1 + y) 
do 
ge (op — ©) u + 3y H;(— w cos ©, 1 — x) 
div ; 
A GE) Hu Sin ©, 4 + 0 cos ow, 1) 
dx 
TS (@o — ©) y 
PRE H;(3u’ Sin ot + 2v sin ot) 
Eu zl< ; 1 A 215 
G Ÿ Il 1 1 (36) 
dy 
a = (o5 — ©) & 
1 out < ! 
+ EY Ha(3u COS @yt— 2u — 3 COS w1 () 
dz : 
re 9+ H,(x sin ©, t + y cos w,t) 
du _ 
dE 


— 3% H,(x Sin wo, t + y cos w;t) | 


(wo — 2Q,, fréquence de résonance). 


Le sous-système (36), contrairement à ce qui se 
passait pour des noyaux de spin / —1, ne se 
réduit pas à 3 équations. Mais on obtient une telle 
réduction si on remplace dans (36) Æ, par sa 
valeur (32), et que l’on fait l’approximation ana- 
logue à celle qui a fait passer de (26) à (27). Le 
sous-système (36) se décompose alors à son tour en 
deux sous-systèmes, dont le premier seul porte sur 
des quantités non nulles : 


mme mo 


N° 8-9 
du 
ANUS @uree0 
do He 
nn Dev. — 0 (37) 
HE, 
dé DA 5 


Les équations (37) sont identiques à celles de 
Bloom, Hahn et Herzog. Il résulte en effet de la 
transformation (33) que 


M, = u cos o,t—vsino,;t+u/ 


1 y# 


ca uz = USMO,t+vcoso;t (38) 


et comme #’ reste nulle au cours du mouvement, 
on voit que les équations (37) sont de la forme (30), 
avec la valeur € — 3 qui correspond à la formule 
de Bloom e. a. (On a pris des valeurs de w, v, et w 
2 fois plus grandes que celles de Bloom e. a. puisque 
on écrit directement les équations pour l’ensemble 
des noyaux et non pour la moitié d’entre eux.) 

Si nous écrivons au lieu de (38) la trans- 
formation : 


M M cos rt Lo sin @1£ 
(39) 


AS? 
2 Que cos &,t 


pl 


5 Qu — = M; sin wit 8 


les quantités qui interviennent dans (37) s’inter- 
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x 1 
Yi, = its 35 


A 
Hs VAS 


Yo,t1 GE 


(9 + 16 cos ©, & + 10 cos 2oit) + Rats 


(— 3 — 2 cos w,t + 5 cos 2w,t) + D 
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prètent simplement, comme des composantes du 
moment quadrupolaire et de l’aimantation : 


dMx 1y# 
ra + (oo e1) (2 À Q%) = 0 


diyi 
RG CE) — (00 — 0) M 


Av 
d o:) 


H,f1y}# :.) 
ee 1 ( x Qu) = 0. 


d 

Ti 
En portant dans ces équations la transfor- 

mation (39), nous obtenons des équations dans le 


système fixe : 


dM; LA le E 
dt ob TE Que 0 


dfiy# ñ 
TE . Qu) — v0 M + 3y H: (5202) = 0 


d 
a | 
Mais, contrairement aux équations de Bloch ou 
aux équations (24-I-bis), valables pour des noyaux 
de spin 1, les équations (41) sont des équations 


approchées, obtenues en négligeant certains termes 
dans les équations (31). 


dy 
& © 


Q:) Y Ho Ê < Qu) = 0. 


5. Résonance quadrupolaire de noyaux de spin 
TI —5/2. — Pour obtenir les équations cherchées, 
on écrit les équations (17) pour Z = 5/2 ; on effec- 
tue ensuite la transformation : 


(4 sin ©, € + 5 sin 2œ,t) 


RER 2 sin wo, t + Sin 2w;t) 
412 V14 

1 
De en (3 
#7 42 Va10 


& COS w,t + cos 2w,; t) 


+ Mie e (4 Sin ©, € + 5 sin 2, t) +- Fi 5 cos w,t + 5 cos 20, t) 


- 2 : ; 1 
Tino tr L Sein 207) + Yan EE (— cos w, t + cos 2o, 1) 


9 VA 


[SA 


y* 24 V3 


ADEME RE A 
on ns VS 


+ Van 3 (4 + cos w, t + 40 cos 2oœjt) + Yi 


: 48 V2 
= UNE 
TA 


10 à: 2 
+ Fu (sin ©, t + 2 sin 2w,t) + Ya 7 6 cos wyt + 2 cos 2w,t) + À PES STAE 


(— 3 — 2 cos w, t + 5 cos 2ow,t) + FER 


: (— 2 sin w,t + sin 2o,t) É M 


OX 4 à : 
Æ VS 4 ——> (— 2 sin w,t + sin 2o,t 
5 19 210 L < 1 t) 
3 V3 à 
—— (— sin ©,t + 5 sin 2o;t) 
1 V5 1 1 
v3 
2 V70 
* 


ci Y5,+1 


(sin &, t + 2 sin 2w;t) 


al 

—— (— 2 + cos w, t + cos 20, t) 
5 VI£ : ? 
30 V2 
7V7 


(— cos &w, t + cos 2o, t) 


(5 sin w, t + sin 2œ;t) 


(&2) 
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ue 


1 Ra drrioe 
= Vaso (5.608 So, 1 + 4) 4 ITS Ste o LEE Ve 
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8 V30 …,*. 5 V30 : 
Y 3 — 4 coS ot + COS 2oit) + 1Yo+ — (— 2 sin &w,t + sin 2oœ,t 
fe | 1 1 ar À ; 1t) 


+ ha D 2 + cos w,t + cos 2œ,t) + han (5 sin w, t + sin 2, t) 


+ Fa 37 (10 + 10 cos ©, t + cos 20, 1) 1 


1 
Yo : (9 cos ©, t + 5 cos 3wit) + ha (8 sin © t + 5 sin 3œ,t) 


14 


ta OK 1 : à 
NE er JD re — (— 3 sin ow,t + sin 3o, 1 
4 24 V7 1 1t) 5 "sv si 1 1 i ) 
Æ y a8 : (3 sin ©, t + 5 sin 3œ,t) + Le (cos w;t + 5 cos 3w, t) 
14 V2 : - = 6 
RE EN ER EE ei 
È = + =) \nuEe Sa AE 
4 UM 1 1 5, TE 1 1 
% 45 Ok 5 V2 3 HONTE 
37 = (== cos wyt + cos So,t) +1 — (— sin w,t + 3sin3ot 
2 vs | 1 1 8 29 vai US 12) 
nn Vase 1 (5 cos &w, t + 9 cos Le lt) + ose (5 sin o;t+3sm3e t) 
% 5 V15 ‘ 4 % 5 V5 


ET ie (—— 3 sin w,t + sin Bot) + Ys49 —— (— cos w, t + cos 30, 1) 


9 


f V5 : , 
+ Vars NS (5 sin &, 4 + 3 sin 3w, 4) +- Vas 5 (5 cos &, t + cos 3o, () 
V \ 
7# * L 

(— 1 + cos 3or1) 


2 V2 ASE 


5 


Le AV 3 ND 
D ÉN ÿ Sin 31 L- Fi cos Soit E Da 2 = Sin 3: 1 
* J 
20e à 10 1 
Poe al + cos 3o,/) + APT 3 sin 3, € “|- Youks g li COS 3, + 5) 


Ve La COS MIO RE 26,1 


Soie EMI TE * 
+ L Fat sin 2 wo, 4 + Ÿ:,44 COS 2 wi 


- VE pOur Ur ONE PERS 


écrire les équations transformées, on posera : 


=mAra0 7 VAS j« __ V210 à 1220 : 7 V15 1* V210 7x \ 
1 439 Pre 288 L51 u sNYh 1 TE 108 Lsi + 20 51/ 
ue V210 V2 
= NT Re 1 x 1 b = NyaK* 
216 72 
2 * VA4 % 5 * VAL x 
o) ere TE v = NrA(EX se + K;) 
7 Ne n(y* V405 x \ . v6 * V105 ,* 
RE YA MT ET Yw) mA ra EPA ie Yo + 55 Yao (43) 
À V7 ç# ; 9 # V7 ,x 
* Dnva(is Re) x TR MR Tes — 5% la 
_5 VE WU Rte US TS 
Y — 36 NY (x, Se "50 Ke VA SNA xà,) 
5 be 3 V5 ,« ; * 3 V5 ,* 
a Énen(is E 25 3) a =-N y PÉRrs 
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Le sous-système utile s’écrit avec ces notations : 


du 
dt 


dv 
dé 
dæ 


dé 


du’ 
dt 


dv’ 


= (205 


dæ Bo, 1 


32 


da 
dt 


32 


(On a omis les quatre dernières équations, qui ne 
seront pas utiles.) (©, = 2Q, est la plus petite des 
deux fréquences de résonance, l’autre est 26,) 


» 
me) 
da avec url 


“ NEURE DU 1 
Si on excite la transition | -} =, -|-- 
D 2) 


champ alternatif 


ot, CH, =, = 0 (15) 


et que l’on remplace, comme précédemment, les 
- termes rapidement variables par leur valeur 
moyenne, (44) se décompose à son tour en plusieurs 
sous-systèmes, dont nous n’écrirons que le premier : 


du 

+ (@y — ©) v'— 0 

© — (@9 — ©) u + LME w = (46) 
dw 1 à KT 

(NA EC 

PAT RINS à FRERE 

Hola 


D’autre part, il résulte de la transformation (42) 
que 
My; = u COS ©, t — v sin w,t 
+ uw’ cos 20, t — v’ sin 2o,t + u” 
(47) 


rh : 
Te Qi: = u sin &,t+ 0 cos ot 


pl= 


+ 2u’ sin 26, t + 20 cos 2; t. 


Or u’, v’ et u” restent nulles au cours du mou- 


1 ; . 
—= Te u Sin @,t + 8u’ sin 26, t — Av cos w,t + 8v’ cos 2, t) 
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= — (op — ©) 0 — 8% H,(w sin ©, 1 + x sin 26, 4— y + y coS 2w,) 
= (op — ©) u — 8y H,(w cos ©, t — x cos 20, 4 + à + y sin 20; 1) 
1 1 Me À 
Fe y H,(2u Sin ©, t — y’ sin 20, t + 20 cos ©, t — v’ cos 2w, 1) 
20,) 0 — y H,{(5w’ sin 2o, 1 — 8x sin ©, 1 — 8y COS w; t) 


= (205 — 20,) u — y H,(5w’ cos 20,1 + 8x cos ©, t — 8y sin w,{) 


301) y 39 y H,(8u’ Sin ©, t — 5u Sin 2,1 — 8v’ cos w, t + 5v cos 2w, t — 5b) 
« 1 : 5 
= = (300 — 30) x —— y H,(8u° cos ©, t — Bu cos 20, t + 80’ sin w, t — 5v sin 2w,t + à) 
4 : ù É 
= — (©9 — ©) y’ + 39 y HA(16u” sin ©, t + Ju — 18b cos 26, t — a sin 26,1 — a’ sin w,1) 


Rae à 
— ©) à + y (Gu” cos wi t — Ju — 18b sin 20, € - a cos 20, L — «4° C0S &, 1). 


vement. Nous pouvons donc remplacer la trans- 
formation (47) par (39), et écrire les équations (46) 
sous la forme : 


dM* 
Fa 


DER 
ÿ 6 5 Q%.) (08 0) Ms 


dt 2 
AL, 1y# # 25, 2 
8 2 [-- MES (cr. Oz: v)) O0 (48) 
d AY A % 25 x ] Y EI 6 Yh x.) 
Dei Res Da: (0) TE LR ONU: 
| h & (co. 98 Qz(0) 2 La 


Si on excite la transition (= 9 L 
champ alternatif 
(49) 


HS =tH y cos20;1!, Hy =; = 10 


on obtient de même quatre équations utiles : 


= + (20, —20,) » = 0 

ue (260— 20) u° + 5e sa (59) 
Lee au = 0 

_ +5 Here 


u, v et u” restent nulles au cours du mouvement 


LL. 0. D cos nas ondes 


à 
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donc la transformation (47) peut être remplacée 
par : 
1y# 


Ms = Mo cos 2oit—7 à UE Sin 20, t (51) 


re 


; : 
— Qyr = MX sin 2o,t FT OE cos 2, 1 


et les équations (50) s’écrivent : 


aMy 
dt 


+ 2o 20  G 0%) — 0 


d'Ay%, + 
( Ë Qu) (205 2&;) Mi 


: Guei)] = 0 (52) 


Rs (e-se0)]- CAC) = 0 


6. Interprétation des équations. — Les équa- 
tions (approchées) de la résonance quadrupolaire 
sont de la forme : 


du ; 

Obs Dai 

dv 

d (oo, wo) u + Cy —w (0) (53) 
de el 

dt CNE 


La signification physique des symboles est 
donnée dans le tableau I, ainsi que la valeur de £. 
w', est la fréquence derésonance, «’, lafréquence du 
champ alternatif. 


TABLEAU I 
SPIN TRANSITION u v w E 
me on —5lle—;ouu] « 
> ( 5 D " :) Mi SL Que : LS [eo : Q.(0)] 8 
Le (LÉ LS un Se 


Le signal d’induction est : 
M, = u c0$ ot vSsin oi tt M, = M, = 0: "154 


Il est donc parallèle au champ alternatif [1], 
contrairement à ce qui se passe en résonance 
magnétique. 

L'intégrale du système (53) qui correspond aux 
conditions initiales : w — 0, v = 0, w — w(0) 
(système de noyaux en équilibre thermique à 
l'instant initial), est : 


LA y H 2 2 
te 


u = #(0) pe (1 — cos opt) 
Ë F ” 
os æ (0) 2 sin Op L (55) 
Op 
Ha) , ts 
de : 1) cos opt + (oÿ — wi)? 
w — w(0) 5 
op 
où 
, H 
oi — (@, — wi)? + (5 | È (56) 
En particulier, à la résonance, 
My = — vsin 6 t — w(0) VEsin wptsin wÿt. (57) 


Calculons l'amplitude du signal d’induction (57) 
à la résonance. Fo est proportionnel à @..(0) 
(tableau I), et Q.(0) est donné en fonction du 
gradient a ‘champ qpar 


Ql0) = 8 q (58) 


où £ est la susceptibilité quadrupolaire. Celle-ci se 
calcule par la même méthode que la susceptibilité 
paramagnétique ; on obtient : 


Br +3)U+01(1—0). 2. 


Nous exprimerons l'amplitude du signal d’induc- 
tion à la résonance en calculant son rapport à 
l'amplitude du signal d’induction en résonance 
magnétique à la même fréquence. Le résultat est 
donné par le tableau IT. 

Les résultats obtenus permettent aussi de 
trouver l’effet produit par une impulsion de champ 
alternatif (méthode d’écho). Pour calculer les popu- 
lations des niveaux il suffit de connaître Q., (et 
aussi, pour Z — 5/2, Y,,). En particulier, à la réso- 
nance, On voit qu'une impulsion de durée r/«, 
(«impulsion de 1800 ») change le signe de w ; ce qui 
permet, en se reportant au tableau I, de calculer la 
valeur de Q., après cette impulsion. 
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TABLEAU II 


TRANSITION 


RAPPORT DE L’AMPLITUDE DU SIGNAL D’'INDUCTION 
A CELLE DU SIGNAL D’'INDUCTION 


EN RÉSONANCE MAGNÉTIQUE A LA MÊME FRÉQUENCE 


SPIN 
4 
3 
2 
5 (25 +7) 
2 D 2 
5 (+5 à 
2 PRET 


7. Conelusion. — On peut se demander si la 
méthode par laquelle nous avons obtenu les équa- 
tions de la résonance quadrupolaire n’est pas trop 
compliquée, et s’il n’aurait pas suffi de partir des 
équations de Bloom, Hahn et Herzog et d’inter- 
préter physiquement les quantités qui y figurent. 
On aurait sans doute pu obtenir ainsi le tableau I, 
du moins pour les valeurs demi-entières du spin. 
Mais la méthode que nous avons utilisée permet 
d'obtenir en plus les résultats suivants : 

1) Les équations (53) sont des équations appro- 
chées, valables lorsque y /1, est suffisamment petit. 
Si on tient compte de l’ensemble des termes, on 
n’est pas seulement amené à rétablir dans les coef- 
ficients des équations (53) les termes rapidement 
variables. On doit aussi (du moins pour 1 > 1) 
introduire de nouvelles grandeurs et de nouvelles 
équations. C’est ce qui explique que les équa- 
tions (53) ne permettent pas de calculer leffet 


‘ 

2 
& = 0,35 
-. = 0,16 


Bloch-Siegert lorsqu'il y a plusieurs fréquences de 
résonance (pour { = 5/2)[13]. 

2) La relaxation, elle aussi, peut avoir pour effet 
d'introduire des termes nouveaux et de nouvelles 
équations. La relaxation spin-milieu, par exemple, 
pour des noyaux de spin 5/2, se traduit [5] par un 
terme dans le second membre de la troisième équa- 
tion (46) qui est une combinaison de w et w’ ; on 
ne peut donc plus éliminer la quatrième équation. 

3) Enfin, les équations obtenues permettent 
l’étude de l'effet Zeeman en résonance quadru- 
polaire. L'existence d’un champ permanent rompt 
les symétries qui avaient permis de réduire le sys- 
tème initial à trois équations ; le sous-système qui 
contient les quantités Y,,, Y,4, ete, fournit lui 
aussi des équations utiles, et celles-ci deviennent 
par conséquent beaucoup plus nombreuses. 


Manuscrit reçu le 7 mai 1958. 


BIBLIOGRAPHIE 


[1] KaszwLer (A.), Experientia, 1992, 8, 1. 

[2] Broom (M.), Haun (E.) et HerzoG (B.), Phys. Rev. 
1955, 97, 1699. 

13] Lurcar (E.), Colloque Ampère 1957, Archives des 
sciences, 1957, 10, fasc. spécial. 

[4] Lurcar (F.), C. R. Acad. Sc., 1954, 238, 1386. 

[5] Lurcar (F.), Thèse, Paris, 1956. 

[6] Waneswess (R.) et Bzocu (F.), Phys. Rev., 1953, 89, 
728 


[7] ROSENFELD (L.), Théorie des électrons, Hermann, 1951, 
p. 22 


[8] Racau (G.), Phys. Reo., 1942, 62, 437. 
[9] Ezxiorr (J.) et Lane (A.), Hand. Physik, 1957, 89, 
390 sq. 

[10] BLocn (K.), Phys. Rev.,1946, 70,460. 

[11] AvanrT (Y.), Théorie de la résonance quadrupolaire, 
in : P. GRIVET e. a., La résonance paramagnétique 
nucléaire, G. N. R.$., 1955. 

[12] Bzocx (F.) et SIEGERT (A.), Phys. Rev., 1940, 57, 522. 

[13] Lurcar (F.), J. Physique Rad. (à paraître). 


LE JOURNAL DE PHYSIQUE ET LE RADIUM 


TOME 19, AOUT-SEPTEMBRE 1958, PAGE 724. 


LA FUSION INDUSTRIELLE EST-ELLE POUR DEMAIN ? 


Par G. VENDRYES, 


Centre d’études nucléaires du C. E; AÀ., Saclay. 


Résumé. — Cet article reproduit une conférence donnée à Paris, le 17 avrik 1958, dans le cadre 
de l’exposition de la Société française de Physique. Il expose les principes fondamentaux sur 
lesquels repose la possibilité de produire de l’énergie par fusion et développe plus en détail les 
perspectives offertes par l’utilisation du phénomène de striction. 


Abstract. — This paper reproduces a lecture given in Paris, on April 17th 1958, during the 


French Physical Society Exhibition. 
the possible use of fusion for producing power. 
discussed in some detail. 


Depuis quelques années, un intérêt considérable 
est accordé aux perspectives offertes par l’énergie 
thermonucléaire. Si cette conquête s'impose à nous 
pour d’impérieuses raisons pratiques, elle constitue 
en même temps une passionnante aventure scienti- 
fique et ouvre un nouveau chapitre de la physique 
que M. Perrin qualifiait récemment d’astro- 
physique expérimentale. Il s’agit bien, en effet, de 
reconstituer sur un petit coin de notre terre des 
conditions analogues à celles qui règnent dans les 
étoiles et d'adapter à notre usage les sources où 
elles puisent l'énergie qu’elles rayonnent. 


Les possiblités du deutérium. -_- Dans cette 
entreprise ambitieuse, le seul avantage dont nous 
disposions est d’avoir, dans une certaine mesure, le 
choix des armes. Pour toute une série de raisons, 
notre combustible sera l'hydrogène ou plus préci- 
sément un isotope rare de l'hydrogène, le deutérium 
Les réactions nucléaires auxquelles nous ferons 
‘ appel sont indiquées sur le tableau I ; elles sont 
connues depuis un quart de siècle. 


Therein are developed the basic principles which underlie 
The approach making use of the pinch effect is 


Cette série de réactions représente finalement la 
formation d’un noyau d’hélium 4 par fusion de 
deux noyaux de deutérium ; cette construction 
s’opère en plusieurs temps et donne naissance à des 
noyaux intermédiaires. 

Dans cette filiation, la réaction fondamentale est 
la première, qui présente l’avantage d’utiliser uni- 
quement un élément stable, existant dans la nature, 
Elle peut s’opérer selon deux voies différentes, dont : 
les probabilités sont comparables. 

Les deux autres réactions, qui utilisent les sous- 
produits dela première — d’une façon grossiè- 
rement analogue à ce qui se passe dans une pile à 
fission, où l’on brûle le plutonium après l’avoir 
formé — sont très importantes aussi. 

La réaction d + & doit son intérêt à sa section 
efficace exceptionnellement élevée, qui la rend plus 
aisée à produire que la réaction du deutérium sur 
lui-même. Elle peut donc, dans une première 
phase, faciliter le démarrage d’une production 
d'énergie de fusion à partir de tritium produit 
artificiellement. De toute façon, lorsqu'on saura 


TABLEAU I 


RÉACTIONS DE FUSION 


sous forme de neutrons 


2 à 
g Sous forme de particules chargées 


& 
& sous forme de neutrons 


+ sous forme de particules chargées 


Roue 
d+d-14 M É LV > 3,6 MeV dont { 
| 2,4 MeV 0,8 MeV | 
d+t— n + 1He °— 17,7 MeV dont 
14,1 MeV 3,6 MeV 
+ He 


d+He#— p 
14 


,1 MeV 3,6 MeV 


—> 18,3 MeV, entièrement sous forme de particules chargées 


| 


4 
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utiliser la réaction d + d, on brülera le tritium au 
fur et à mesure de sa production et bénéficiera 
d’un appoint énergétique très important. 

Quant à la dernière réaction d + °He, elle pré- 
sente l’inconvénient d’une section efficace plus 
faible, due à la double charge du noyau d’hélium, 


-mais par contre a le gros intérêt de dégager une 


énergie très élevée, et ce entièrement sous forme de 
particules chargées qui peuvent donc être conser- 
vées dans le milieu réactionnel et contribuer à 
l'entretien de son énergie interne. 

La consommation intégrale d’un kilogramme de 
deutérium par la réaction d + d seule dégagerait 
24 millions de kWh, production qui se trouve être 
équivalente à ce que donnerait la consommation 
intégrale d’un kilogramme d’uranium ou de plu- 
tonium par fission, et qui est environ 3 millions de 
fois supérieure à l’énergie fournie par la combustion 
complète d’un kilogramme de bon charbon. 

Si nous tenons compte de la réutilisation du 
tritium formé, ce qui est légitime, un kilogramme 
de deutérium peut donner 67 millions de kWh et 
une centaine si le He est également consommé. 

On peut difficilement rêver d’une meilleure 
source d’énergie. Si nous songeons à l'immense 
réservoir que constituent les océans, nous trouvons 
là de quoi assurer les besoins de la planète pendant 
des millions d’années. 


Conditions à satisfaire. — Actuellement, la seule 
possibilité qui soit sérieusement envisagée pour 
libérer à grande échelle cette énergie de fusion 
consiste à porter le gaz deutérium à une tempé- 
rature très élevée, de telle sorte que l’énergie 
d’agitation thermique des noyaux devienne d’un 
ordre de grandeur suffisant pour leur permettre, 
lors de leurs rencontres, de surmonter les forces de 
répulsion - électrostatique et de donner lieu aux 
réactions nucléaires attendues. 

Divers procédés, certains fort subtils, ont été 
proposés qui auraient éventuellement permis d’uti- 
liser industriellement les réactions de fusion à la 
température ordinaire. Je mentionnerai seulement 
l’idée, émise voilà quelque temps, de provoquer des 
réactions en chaîne de fusion en utilisant comme 
catalyseurs, comme véhicules de la réaction, des 
mésons {& négatifs de même que les neutrons pro- 
pagent les réactions de fission dans une pile à 
uranium. Les espoirs un moment éveillés par cette 
fusion froide se sont avérés non fondés, non pas 
tant pour des raisons de principe, que pour des 
questions d'ordre de grandeur. D’autres formules, 
reposant par exemple sur l'emploi d’accélérateurs, 
se heurtent à des objections semblables. 

Sans qu’on puisse totalement exclure l’existence 
d’un raccourci, il est clair actuellement que le 
chemin aboutissant à la domestication de l’énergie 
de fusion nous emmène d’abord vers la quête des 
températures très élevées. 

Pour fixer les idées, nous allons partir d’une 
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certaine quantité de deutérium pur à la tempé- 
rature ordinaire et à basse pression, égale par 
exemple au dix millième de Îa pression atmo- 
sphérique — ce qui correspond à un nombre dé 
noyaux de deutérium de l’ordre de 5.10% par en° 
— et nous supposons que nous chauffons ce gaz à 
volume constant sans nous soucier pour le moment 
des problèmes posés par cette opération. Quelle sera 
la succession des phénomènes observés ? 

Vers 5 000 °C, les molécules vont être disloquées 
en atomes, puis bien avant que le cap des 
100 000 ©C soit atteint, les atomes de deutérium 
vont être dissociés en leurs constituants, noyaux et 
électrons, le gaz complètement ionisé constitue 
alors ce qu’il est convenu d’appeler un plasma. 
Nous continuons à chauffer. C’est seulement une 
fois que le million de degrés est dépassé que des 
appareils très sensibles peuvent permettre de 
détecter les rares neutrons attestant que des réac- 
tions de fusion commencent à se produire au cours 
de collisions exceptionnelles entre noyaux de deu- 
térium. Lorsque la température augmente, le taux 
de réactions thermonucléaires croît très vite, la 
puissance dégagée aussi. À 100 millions de degrés, 
elle atteint 10 watts par cm, tandis que la 
pression exercée par le plasma dépasse 100 atmo- 
sphères. Nous ne sommes toutefois pas encore 
arrivés au but car, dans de telles conditions, l’éner- 
gie qui s'échappe du plasma sous forme de rayon- 
nement électromagnétique, dû surtout au freinage 
des électrons dans le champ des noyaux, excède 
encore l’énergie produite par les réactions de fusion. 
Il s’ensuit qu’à moins d’un apport continu d'énergie 
extérieure, la température du mieu ne peut que 
décroitre et le dégagement d’énergie par fusion 
s’éteindre progressivement. 

C’est seulement lorsque la température atteint 
400 millions de degrés que l’énergie produite par 
fusion devient supérieure à celle qui est rayonnée ; 
à condition qu'aucune autre cause de pertes n'inter- 
vienne, le dégagement d’énergie peut alors s’entre- 
tenir de lui-même, comme dans une combustion 
chimique ordinaire. Il est même vraisemblable 
qu’un effort supplémentaire sera nécessaire pour 
compenser d’autres causes de pertes, par conduc- 
tion notamment, et que la température de fonc- 
tionnement devra se tenir plutôt aux environs du 
milliard de degrés, auquel cas la contribution de la 
réaction d + *He deviendra essentielle. Dans ces der- 
nières conditions, la pression de notre plasma 
atteindra 1 000 atmosphères et la puissance dégagée 
500 watts par cm. 

Une analyse de ce genre montre que les princi- 
paux paramètres physiques de fonctionnement 
d’un milieu thermonucléaire au deutérium prati- 
quement utilisable ne peuvent s’écarter d’une 
gamme de valeurs relativement étroite. Comme 
nous venons de le voir, la température à atteindre 
doit dépasser un seuil critique de quelques centaines 
de millions de degrés. Quant à la densité du milieu, 
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nous ne pouvons la prendre très différente de celle 
du précédent exemple (soit 1015 à 1016 noyaux 
par em?) sans tomber dans des ordres de grandeur 
inadmissibles pour la puissance produite ou la 
pression à contrebalancer. 

Dans ces conditions, le temps nécessaire à la 
consommation d’une fraction notable du deutérium 
présent est de l’ordre d’une dizaine de secondes. Il 
n’est pas inutile de signaler que la trajectoire par- 
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F1G. 1.— La courbe A représente la puissance produite 
par fusion en fonction de la température dans un plasma 
constitué uniquement de deutérium et contenant 5.1015 
noyaux par Cm?. 

La courbe B représente la même grandeur dans le cas 
d’un plasma contenant par em? 2,5.1015 noyaux de 
deutérium et 2,5.10'5 noyaux de tritium. 

Enfin la courbe C représente la puissance qui s’échappe 
sous forme de rayonnement électromagnétique d’un 
plasma hydrogéné complètement ionisé contenant 5.1015 
noyaux par Cm. 

es ordonnées de ces courbes sont toutes trois propor- 
tionnelles au carré de la densité des noyaux dans le 
plasma ; aussi les températures critiques pour lesquelles 
les courbes A et B coupent la courbe C sont-elles indé- 
pendantes de cette densité. 
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courue pendant ce temps dans le plasma par un 
noyau individuel a une longueur comparable à 
celle de la circonférence terrestre. 

Dans le cas où nous utiliserions comme combus- 
tible un mélange à parts égales de deutérium et de 
tritium, les conditions seraient un peu moins 
sévères. La température critique est alors seule- 
ment de 50 millions de degrés. A titre indicatif, un 
plasma contenant 5.101# noyaux de deutérium 
et 5.107 noyaux de tritium par cm°, porté à 
400 millions de degrés, exercerait une pression de 
100 atmosphères, dégagerait une puissance de 
500 watts par em et se consumerait en un temps 
de l’ordre de la seconde. 

Pour atteindre les objectifs qui viennent d’être 
précisés, deux problèmes fondamentaux sont à 
résoudre. 11 faut trouver les moyens de porter le 
plasma aux hautes températures requises, il faut 
aussi le confiner, c’est-à-dire l’isoler des parois du 
récipient qui l’enferme. Ces deux problèmes sont 
intimement liés, car on doit les résoudre simul- 
tanément, et en faisant appel aux mêmes pro- 
priétés physiques du plasma. Pour la clarté de la 
présentation, il est bon de les évoquer succes- 
sivement. ; 


Le problème du chauffage. — Pour porter à un 


milliard de degrés un litre de deutérium se trou- 
vant initialement à la température ordinaire, sous 
une pression égale au dix milième de la pression 
atmosphérique, il faut lui fournir une énergie 
équivalente à la quantité de chaleur qui est néces- 
saire pour chaufter de 50 0C un litre d’eau. Une 
telle énergie est sans doute modeste, en valeur 
absolue ; elle représente néanmoins une fraction 
non négligeable (quelques pour cent) de l’énergie 
maximum qu’il serait possible d'extraire par réac- 
tions de fusion du plasma chaud, dans le cas 
extrême où l’on parviendrait à consommer ainsi 
l'intégralité du deutérium et de ses dérivés, tré 
et hélium 3. 

En réalité, il nous faudra mettre en œuvre une 


énergie nettement plus importante afin de com-. 


penser, pendant tout le temps que durera le chauf- 
fage, les pertes — principalement les pertes par 
rayonnement —. Nous avons, en quelque sorte, à 
remplir à ras bord un tonneau qui fuit. On conçoit 
que, pour y parvenir, il faille opérer aussi vite que 
possible. Dans le cas de notre plasma, celà signifie 
que nous devrons utiliser pour le chauffer une 
puissance très grande. Pour traverser la région 
située en dessous de 100 000 °C, où le gaz incom- 
plètement ionisé rayonne énormément et présente 
une capacité calorifique élevée, la puissance à 
mettre en œuvre devra même atteindre pendant de 
brefs instant, des valeurs considérables. 

En soi, le problème du chauffage ne devrait pas 
soulever de difficultés d'ordre fondamental. On sait 
que la chaleur est la forme la moins noble de 


. 
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l'énergie, celle en laquelle toutes les autres ont une 
tendance spontanée à se transformer, la difficulté à 
laquelle on se heurte dans la plupart des cas, étant 
justement d'éviter ou de limiter ces transfor- 
mations. [ci nous nous trouvons au a contraire avoir 
à les favoriser. 


L'un des principaux parmi les processus irréver- 
sibles est l'effet Joule et c’est effectivement lui 
qu’on met à profit dans les dispositifs actuellement 
étudiés, en faisant circuler dans le plasma un 
courant électrique intense. Toutefois, la résistivité 
d’un plasma hydrogéné complètement ionisé 
_décroît selon une loi en T— *”? lorsque la température 
augmente, devenant inférieure à celle du cuivre au 
delà de 30 millions de degrés. De plus, ce sont 
essentiellement les électrons qui transportent le 
courant et qui transfèrent ensuite aux ions, par 
des chocs, l'énergie que le champ électrique leur a 
communiquée. Or ce sont les ions que nous désirons 
chauffer ; nous sommes donc là dans des conditions 
relativement défavorables et qui le deviennent de 
plus en plus au fur et à mesure que la température 
s'élève et que le transfert d'énergie des électrons 
aux ions se ralentit sans cesse. En conclusion, 
l’effet Joule sera surtout efficace pendant les pre- 
mières phases du chauffage, très intéressant notam- 
ment pour passer le cap des 100 000 o0C. 


Un autre procédé bien connu pour élever la 
température d’un gaz consiste à le comprimer. Il 
est également utilisable pour chauffer un plasma, 
la compression étant alors provoquée par les efforts 
électro-dynamiques s’exerçant sur les courants 
présents dans le plasma. Si ces compressions sont 
assez brutales, elles donnent naissance à de véri- 
tables ondes de choc d’autant plus efficaces pour le 
but que nous poursuivons qu’elles augmentent 
lP’énergie cinétique des ions de préférence à celle des 
électrons. L'application répétée de compressions 
irréversibles, en cumulant leurs effets, offre ainsi un 
puissant moyen de chauffage. 


Il en est d’autres encore susceptibles de prendre 
le relais de l’effet Joule aux températures élevées, 
et qui reposent, sous une forme ou une autre, sur la 
transformation d'énergie électromagnétique en 
chaleur dans le plasma. 


Il est toutefois un dernier procédé que je voudrais 
mentionner car il est très différent de tous les 
autres et peut éventuellement apporter une solu- 
tion nouvelle au problème du chauffage. Il con- 
siste à alimenter directement le milieu réactionnel 
en noyaux et électrons accélérés sous quelques 
dizaines de milliers de volts, issus de sources 
d’ions à forte intensité. Le plasma chaud se cons- 
titue par simple brassage de ces particules, au 
cours duquel leur énergie moyenne demeure 
inchangée. La formule est élégante mais les inten- 
sités requises dépassent ce que peuvent donner les 
sources d’ions actuelles. 
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Le problème du confinement. — L'autre pro- 
blème fondamental que nous avons à résoudre est 
celui du confinement. Comment inventer un réci- 
pient qui contienne le plasma à des pressions et à 
des températures aussi élevées ? Il est bien évident 
qu'aucune paroi matérielle ne peut supporter des 
températures dépassant quelques centaines de 
degrés. Il faut donc trouver les moyens d’établir un 
gradient de températures considérable entre le 
cœur du plasma et la surface interne du récipient 
qui le renferme. 

Une formule possible, sur le papier, consisterait 
à construire une grande sphère métallique, à la 
remplir entièrement de gaz, à en porter le centre 
aux températures nécessaires à la fusion et à y 
établir un régime tel que l’énergie dégagée dans la 
région centrale soit évacuée principalement par 
conduction à travers des couches gazeuses de plus 
en plus froides jusqu’à la paroi maintenue à la 
température ordinaire. Le calcul d’un tel engin 
méritait d’être effectué, mais ses conclusions en 
condamnent la réalisation pratique, même dans le 
cas estimé le plus favorable où l’on envisage d’uti- 
liser comme combustible un mélange de deutérium- 
tritium sous une pression de 1 000 atmosphères. La 
conductibilité thermique extrêmement grande du 
plasma chaud obligerait en effet à donner à cette 
étoile miniature un rayon de { kilomètre tandis que 
la puissance produite dépasserait 1017 watts. 

Il faut donc admettre que dans tout appareil 
réalisable à l’échelle terrestre, le plasma sera prati- 
quement isotherme, et devra être strictement 
maintenu hors de tout contact avec les parois du 
récipient. [l y a lieu de préciser que cette exigence 
vise beaucoup plus à protéger le plasma que la 
paroi. En cas de contact entre eux, la paroi 
s’échaufferait modérément, alors que le refroidis- 
sement du plasma serait considérable, à cause des 
pertes par conduction et également par suite de la 
pollution du gaz par des particules de numéro 
atomique élevé arrachées à la paroi. 

Pour réaliser l’isolement total du plasma dans 
un récipient de dimensions raisonnables, il semble 
bien que nous ne disposions que d’un seul moyen. 
Il consiste à utiliser le fait que le gaz est complè- 
tement ionisé et que ses particules constitutives 
sont donc sensibles à l’action des forces électro- 
magnétiques. Dans un champ magnétique en parti- 
culier, elles décrivent des trajectoires enroulées en 
hélice le long des lignes de force et peuvent ainsi 
être maintenues dans une région finie de l’espace. 
On cherchera donc à enfermer le plasma au moyen 
d’une barrière immatérielle en remplissant l’espace 
qui l’entoure avec du champ magnétique. Suppo- 
sons le cas idéal représenté sur la figure 2. où nous 
aurions du champ à l’extérieur et pas de champ à 
l'intérieur du plasma ; les particules emprisonnées 
viennent rebondir à la limite du champ magnétique. 
Les déplacements de ces charges positives et néga- 
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tives, au moment où elles sont courbées par le 
champ à la frontière du plasma, donnent naissance 
à un courant électrique superficiel. On peut égale- 
ment, d’un point de vue macroscopique, se repré- 


plasma FA 
PER ES A E 


Fic. 2. — Schéma du confinement d’un plasma complète- 
ment ionisé par un champ magnétique qui l'entoure ; les 
flèches représentent le courant électrique superficiel dû 
aux particules courbées à la lisière du champ magnétique. 


senter que le confinement du plasma est dû à la 
force exercée sur ce courant par le champ extérieur, 
et dirigée vers l’intérieur. La pression du gaz 
enfermé est ainsi contrebalancée par une pression 
magnétique antagoniste, proportionnelle au carré 
du champ. Pour équilibrer une pression de 
300 atmosphères, un champ de 100 000 Gauss est 
nécessaire. Ce sont évidemment là des valeurs con- 
sidérables, mais qui ne semblent pas inaccessibles. 

Pour réaliser une telle bouteille magnétique, 
diverses configurations spatiales de champs magné- 
tiques peuvent être imaginées, dues à des courants 
électriques circulant dans le plasma ou dans des 
conducteurs extérieurs. 

Finalement, on constate que les deux objectifs 


poursuivis, chauffage et confinement, peuvent au 


moins en principe être atteints en faisant appel 
à des méthodes de même nature. 


. Le phénomène de striction. — L'un des procédés 
actuellement envisagés utilise ce qu’on appelle le 
« pinch effect » dans les pays anglo-saxons, et que 
nous nommons le phénomène de striction. I] s’agit 
là d’un procédé parmi d’autres, et pas nécessai- 
rement le meilleur bien que de grands efforts lui 
soient aujourd’hui consacrés. Si je m°y arrête, c’est 
qu'il est le seul sur lequel ait été levé le secret qui 
malheureusement entoure encore la plupart des 
travaux sur la fusion thermonucléaire. 

Ce phénomène, connu depuis longtemps, est 
représenté schématiquement sur la figure 3. On le 
reproduit aisément au laboratoire en déchargeant 
une batterie de condensateurs dans une chambre 
remplie de gaz sous basse pression, entre élec- 
trodes. Le courant qui circule dans le gaz, paral- 
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lèlement à l’axe de la chambre, crée un champ 
magnétique azimutal dont leffet, sur le courant 
qui le produit, est une force dirigée vers l’axe. Si le 
courant est assez intense, —- dépassant dañs un 
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F1G. 3. — Représentation schématique d’une décharge 
entre électrodes contractée sur l’axe de la chambre par 
le phénomène de strietion. k 


cas type 200 000 ampères — la force ainsi exercée 
est suffisante pour comprimer le gaz et le rassem- 
bler dans une colonne axiale, isolée des parois de la 
chambre. Nous avons donc là un moyen de chauffer 
le plasma, tant par l'effet de compression que par 
l’eflet Joule et de l’enfermer dans un corset de 
lignes de forces magnétiques circulaires, perpen- 
diculairement à l’axe. 

Les premières déclarations relatives à l’appli- 
cation de ces décharges à forte intensité entre 
électrodes aux recherches sur la fusion furent faites 
en mai 1956 par l’Académicien soviétique Kur- 
chatov. Elles eurent un grand retentissement et 
furent suivies par la publication de travaux ana- 
logues entrepris dans divers pays, notamment aux 
U.S. À., en France, en Angleterre, en Suède, etc. 

Les températures atteintes au cours de telles 
décharges dans du deutérium ne semblent pas avoir 
dépassé notablement 100 000 °C, pendant des 
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durées de l’ordre de la microseconde. Quant aux 
neutrons qui sont observés au moment où la 
colonne est la plus resserrée, il est maintenant 
établi, à la suite de très concluantes expériences 
américaines, qu’ils ne sont pas d’origine thermo- 
nucléaire mais que les réactions de fusion aux- 
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Fic. 5. — Dispositif de décharges loroïdales installé au 
Centre d'Études Nucléaires de Fontenay-aux-Roses. Le 
tore en cours d’essai est en verre et comporte une enve- 
loppe extérieure en cuivre sur laquelle sont bobinés les 
circuits inducteur et stabilisateur. 

On distingue des enroulements destinés à la mesure de 
l’intensité du courant dans la décharge. Sur la droite, 
un tore de démonstration en verre, pourvu de huit 
-noyaux magnétiques. 


comme moyen d'investigation, elles ne semblent à 
aucun titre pouvoir constituer l’amorce d’un réac- 
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quelles ils doivent leur origine sont provoquées par 
des deutons accélérés parallèlement à l’axe de la 
décharge dans des champs électriques intenses qui 
se forment localement et de façon fugitive, dans les 
replis de la colonne du plasma. 

Si de telles décharges sont fort intéressantes 


circuit magnétique 


; à a de courent 
/ elecirIque 


Lignes de Champ mégnétique 


4. — Principe de la décharge annulaire, sans électrodes. 


teur à fusion, la présence des électrodes empêchant 
l'isolement complet du plasma. 

Aussi depuis plusieurs années, et d’abord en 
Angleterre et aux Etats-Unis, des recherches ont- 
elles été entreprises vers la réalisation de décharges 
gazeuses annulaires, refermées sur elles-mêmes, 
sans électrodes. Les nombreux dispositifs construits 
à cette fin fonctionnent comme de véritables trans- 
formateurs dont le secondaire est fait d’une spire 
unique : l’anneau de plasma maintenant complè- 
tement isolé dans l’espace (cf. fig. 4). 

Les travaux entrepris sur ces décharges droites 
ou toroïdales mirent très vite en évidence un phé- 
nomène très désagréable. Si aucun dispositif n’est 
prévu pour y parer, la colonne de plasma se tord, 
se disloque, se fragmente en un temps très court 
qui peut être de l’ordre de la microseconde et vient 
frapper les parois, mettant ainsi fin à l’entreprise. 

Parmi les divers procédés proposés pour remédier 
à cet état de chose, les deux principaux sont 
l'emploi d’un champ magnétique longitudinal créé 
à l’intérieur de la colonne de la décharge par des 
bobinages extérieurs et l'emploi d’une enveloppe 
conductrice qui entoure le tube à décharge, et agit 
grâce aux courants induits qui y sont le siège 
lorsque la colonne de plasma se déforme et se 
déplace à son intérieur. Le champ magnétique 
axial est surtout efficace contre les déformations 
dont l’échelle de longueur, parallèlement à l’axe, 


est petite vis-à-vis du diamètre de la colonne ; 
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l’effet de l'enveloppe conductrice est, par contre, 
d'autant meilleur que cette échelle de longueur est 
plus grande vis-à-vis du diamètre intérieur de 
l'enveloppe. On peut donc espérer, par une combi- 
naison des deux, tendre à une stabilisation com- 
plète de la colonne, la question étant plutôt de 
savoir si les conditions assurant cette stabilisation 
peuvent être maintenues pendant les durées de 
plusiéurs secondes requises pour le fonctionnement 
d’un engin thermonucléaire à fusion. 

La manière pratique de lutter victorieusement 
contre l’instabilité foncière de la colonne de la 
décharge, tout en assurant le chauffage du plasma, 
constitue l’objectif principal des appareils récem- 


Fi6. 6. — Dispositif de décharges toroïdales installé par le 
C. E. À. au Centre d’Essais de l’Électricité de France à 
Fontenay-aux-Roses. Le tore essayé ici est en aluminium; 
on distingue sur la photographie les éléments du cireuit 
magnétique améliorant le couplage entre le circuit induc- 
teur et la décharge, les bobines créant le champ magné- 
tique stabilisateur et divers appareils (poste HF pour la 
préionisation du gaz, spectrographe optique, appareil de 
cinématographie ultra-rapide, etc..). 


ments décrits tels que le Perhapsatron S3 de Los 
Alamos ou la machine Zeta de Harwell. Les 
figures 5 et 6 reproduisent la photographie de deux 
montages de décharges toroïdales, actuellement 
essayés par le Commissariat à l'Énergie Atomique. 
Tous ces appareils sont essentiellement des outils 
expérimentaux et sont à coups sûr aussi éloignés des 
réacteurs thermonucléaires de l’avenir que les 
premières machines à vapeur le sont des turbines 
modernes. Il n’est pas cependant exclu qu’ils con- 
tiennent dès maintenant, même dans un état 
embryonnaire, l’essentiel des dispositifs qui per- 
mettront plus tard de faire fonctionner des engins 
industriels. 


Visions d’avenir. — Toute extrapolation ne peut 
aujourd’hui comporter qu’une grande part de fan- 
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taisie. Si, cependant, nous essayons d'imaginer ce 
que pourrait être un réacteur thermonucléaire uti- 
hsant le principe de la striction —— en admettant 
que l'instabilité de ce phénomène ne constitue pas 
un obstacle insurmontable— qu’entrevoyons nous ? 

Nous imaginons un tore de grande dimensions, 
imposées plus encore par des considérations techno- 
logiques que pour des raisons à caractère fonda- 
mental. Son petit diamètre serait au minimum de 
l'ordre du mètre, son grand diamètre supérieur 
à 10 mètres. Le récipient pourrait être constitué 
d’une épaisse enveloppe en aluminium dont on 
aurait rendu la paroi interne isolante en formant 
par exemple une couche superficielle d’alumine. 
Autour seraient disposés les bobinages créant le 
champlongitudinalstabilisateur et le circuit magné- 
tique assurant le couplage entre les circuits induc- 
teurs et la décharge. Nous pouvons par exemple 
envisager un fonctionnement cyclique, en renou- 
velant à chaque cycle la charge de deutérium intro- 
duit dans l’appareil. Dans un premier temps, 
pendant lequel le gaz est ionisé et porté à plus de 
500 millions de degrés, une puissance électrique 
de 1 000 MW est nécessaire, pendant une durée de 
l’ordre de la seconde, pour chauffer un plasma et 
remplir l’espace avec les champs magnétiques 
assurant le confinement stable. Le chauffage s’effec- 
tuerait d’abord par effet Joule, puis au moyen 
d’une série de compressions exercées sur le plasma 
en pulsant le courant inducteur ou le courant 


créant le champ stabilisateur. Pendant le temps 


suivant, qui doit durer plus d’une dizaine de 
secondes, un courant de l’ordre de 10 millions 
d’ampères se maintient dans le plasma très con- 
ducteur, tandis que les réactions de fusion y 
dégagent une puissance de l’ordre de 300 MW. 

La possibilité a été mentionnée à diverses 
reprises qu'un réacteur à fusion permette la pro- 
duction directe d'électricité, sans passer par l’inter- 
médiaire de la chaleur. Le schéma auquel on songe 
à ce propos est le suivant : lorsque des réactions de 
fusion se développent dans le plasma, il s’'échauffe, 
sa pression augmente, il tend donc à distendre/le 
corset de champ magnétique qui l’enserre et cette 
action contre les forces électro-dynamiques va se 
traduire par l’apparition dans des bobinages exté- 
rieurs de courants induits. Il faut cependant bien 
insister sur le caractère nécessairement très vague 
de tels schémas, qui ne constituent à l’heure 
actuelle qu’une perspective séduisante. Dans tous 
les cas il est bien clair qu’une part très notable de 


l'énergie produite sera inéluctablement convertie - 


dans les régions extérieures de l’appareilen chalenr 
à basse température, soit par conduction du plasma 
à la paroi, soit par l'intermédiaire du rayonnement 
électromagnétique et des neutrons s’échappant du 
milieu réactionnel. 

À puissance produite égale, les neutrons issus 
d'un réacteur à fusion au deutérium seront 
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plusieurs dizaines de fois plus abondants que ceux 
qui sont émis par le cœur d’une pile atomique à 
fission de luranium. Ils devront eux aussi être 
absorbés dans un épais blindage extérieur, où ils 
pourront être utilisés pour régénérer par exemple 
du tritium à partir de lithium 6. L'existence d’un 
rayonnement aussi intense posera à coup sûr de 
graves problèmes pour le fonctionnement et l’entre- 
tien de l’appareillage entourant l'enceinte où 
s’opèrent les réactions de fusion. S'il est vrai que les 
produits directs de ces réactions sont incompa- 
rablement moins actifs que les produits de fission 
de l’uranium, on constate cependant que les réac- 
teurs à fusion ne seront pas exempts des problèmes 
techniques posés par la radio-activité. 

On peut, par contre, être tout à fait tranquille en 


ce qui concerne les risques présentés par le fonc- 


tionnement d’un réacteur thermonucléaire. Il 
n'existe en effet aucune façon qu’un tel réacteur 
puisse s’emballer dangereusement ou donner lieu à 
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quoi que ce soit ressemblant à une explosion 
nucléaire. 

Les considérations qui précèdent ont, je 
pense, suffisamment mis en évidence les énormes 
difficultés que soulève la mise au point d’un réac- 
teur thermonucléaire viable et la distance qu’il 
nous reste à parcourir pour y parvenir. Aussi à la 
question : «La fusion industrielle est-elle pour 
demain ?» devons-nous répondre : Non, en tout cas 
pas pour demain matin. Cependant, le fait 
qu'aucune impossibilité fondamentale n'ait été 
reconnue et l’ampleur des moyens qui seront sans 
nul doute consacrés à une entreprise dont l'enjeu 
est aussi important doivent nous inciter à l’opti- 
misme. Aussi voudrais-je en concluant exprimer la 
conviction que la fusion industrielle sera la réalité 
d’après-demain, et si ce ne doit être que la semaine 
prochaine je ne pense pas que nos petits-enfants 
puissent nous en vouloir. 

Manuscrit reçu le 3 juin 1958. 
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CONSIDÉRATIONS 
SUR LES SOURCES D'’'ÉTINCELLES 
DANS LE VIDE E 
POUR + 
ANALYSE SPECTROCHIMIQUE D'ÉMISSION 
DANS L'ULTRAVIOLET LOINTAIN 


Par MM. J. Romaxp et B. Vopar, 
Laboratoire des Hautes Pressions (Bellevue). 


Dans des travaux antérieurs [1], [2] nous avons mis 
en évidence la possibilité d'utiliser à des fins ana- 
lytiques le spectre d’atomes fortement ionisés émis 
dans l’ultraviolet lointain par les étincelles dans le vide. 

Mais l’étincelle dans le vide, sous sa forme classique 
dite «étincelle de Millikan » ne répond absolument pas 
aux conditions de l’analyse spectroscopique, car l’espa- 
cement des électrodes y est réduit à une fraction de 
millimètre et cela occasionne une instabilité prononcée 
de son fonctionnement. 

La source que nous avons particulièrement étudiée 
et utilisée est l’étincelle glissante dans le vide [3], [4], 
qui supprime cet inconvénient par suite de l’abais- 
sement du potentiel disruptif. Mais elle introduit 
comme support de l’étincelle un matériau étranger 
dont les particules participent, au moins partiellement, 
à l’émission, ce qui peut constituer un inconvénient 
majeur dans certains cas, comme par exemple, l’étude 
des possibilités de détection des gaz contenus dans le 
matériau des électrodes. 

La source la plus satisfaisante serait donc une 
étincelle dans le vide, sans support, mais le fonction- 
nement plus commode que l’étincelle classique. 

Pour créer le bombardement initial, nécessaire à 
l’amorçage d’une décharge sous vide, nous avons 
adopté actuellement le procédé qui consiste à utiliser 
une étincelle auxiliaire, procédé qui a l'avantage d’être 
une réalisation simple. — Ce procédé a été utilisé pour 
provoquer l’émission de rayons X, mais on ne savait 
pas s’il pouvait s'appliquer à la réalisation d’une 
source spectroscopique, notamment pour l’ultraviolet 
lointain. Nous avons constaté qu’une telle décharge 
sous vide avec étincelle auxiliaire constituait effecti- 
vement, dans des conditions convenables, une source 
émettant avec intensité des raies d’atomes fortement 
ionisés, caractéristiques du matériau d’une électrode. 

La chambre à étincelles renferme trois électrodes, la 
cathode étant commune au circuit de l’étincelle auxi- 
liare et à celui de l’étincelle principale. L’étincelle 
auxiliaire se produit entre deux électrodes coaxiales 
dont l’écartement est de quelques dixièmes de milli- 
mètres, car cette étincelle fonctionne sans support ; il 
est d’ailleurs aisé de prévoir une disposition permettant 
la formation d’une étincelle glissante, donc l’utilisation 
d’une tension d’amorçage plus basse. L’électrode à 
analyser -(anode du circuit principal) est disposée en 
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face de cette étincelle auxiliaire et a la forme d’un 
crayon métallique dont la forme de l’extrémité n’est 
pas critique. 

La puissance dépensée dans le circuit de l’étin- 
celle auxiliaire peut être relativement faible 
(0,1 &F — 20 KV), par contre il est essentiel d'utiliser 
pour l’étincelle d’analyse une forte capacité (1 ou2uF). 
En effet, pour de faibles-puissances la vaporisation du 
matériau de l’anode est faible et dans ce cas, évidem- 
ment, on observe surtout une forte émission de 
rayons X. De plus, il y a intérêt à favoriser la phase de 
la décharge se produisant dans la vapeur de l’anode, 
relativement à celle qui se produit dans les gaz rési- 
duels, si l’on veut obtenir particulièrement en ce qui 


concerne les gaz, un spectre provenant exclusivement 


du matériau de l’anode. 

Les spectres que nous avons obtenus dans ces condi- 
tions sont analogues à ceux donnés par l’étincelle glis- 
sante et sont caractéristiques du matériau de l’anode 
principale. Toutefois, pour les fortes énergies d’exci- 
tation (fortes capacités) on observe une croissance de 
l'intensité des raies de forte excitation en fonction de la 
capacité, ce dernier résultat étant différent de celui 
observé pour les étincelles glissantes [5]. De plus, le. 
spectre continu qui apparaît pour les fortes excitations 
est beaucoup plus faible que dans le cas des étincelles 
glissantes. Ces deux résultats pourraient s’interpréter 
en admettant, ce qui est vraisemblable, que la densité 
du plasma est plus faible au voisinage de Panode que 
dans l’espace relativement réduit compris entre les 
deux électrodes de l’étincelle glissante. 

Au point de vue pratique, les raies sensibles des 
métalloïdes contenus à faible concentration dans le 
matériau de l’anode principale apparaissent (nos essais 
actuels ont porté sur le soufre, le phosphore, et le 
carbone), et cette source peut être utilisée à la détee- 
tion et au dosage de ces éléments. On observe éga: M 
lement des raies sensibles appartenant au spectre de M 
l’oxygène et de l’azote ; nous avons constaté que l’in- 
tensité de ces raies est bien une fonction de la concen- 
tration de ees gaz dans l’anode et nous étudions actuel- 
lement les possibilités effectives de cette source pour 
le dosage des gaz contenus dans les matériaux solides. ;. è 


Lettre reçue le 25 juin 1958. 
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